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abstract 


The  primary  goal  of  this  thesis  was  the  development  of  an 
experimental  technique,  and  supporting  software,  for  the  ac¬ 
quisition  and  analysis  of  unsteady  velocity  data  generated  by 
an  oscillating  airfoil.  This  research  was  in  support  of  a 
major  investigation  of  the  compressibility  effects  on  dynamic 
sta  1 1  . 

The  experimental  procedure  involved  schlieren  flow  visu¬ 
alization  for  comparison  of  steady  and  unsteady  flow  fields, 
and  for  determination  of  parameters  for  further  study.  Laser 
Doppler  veloc iraetry  was  employed  for  obtaining  velocity  data 
in  the  airfoil  wake.  For  unsteady  data,  the  airfoil  was 
oscillated  in  pitch  about  its  quarter  chord. 

The  data  analysis  produced  wake  profile  plots  represent¬ 
ing  the  flow  field  disturbed  by  the  airfoil.  Results  were 
obtained  for  steady  and  unsteady  conditions. 
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A.  DYNAMIC  STALL 

"he  possible  benefits  and  consequences  inherent  in 
dynamic  stall  have  stimulated  experimental  investigations 
leading  to  significant  progress  in  understanding  this  complex 
p h e n om enon.  Much  of  the  earlier  work  in  this  area  was  car¬ 
ried  cut  by  McCroskey  using  a  two-dimensional  airfoil 
section  oscillating  in  pitch  about  its  quarter-chord.  The 
results  of  this  effort  succeeded  in  clarifying  the  features 
of  dynamic  stall  up  to  a  Mach  number  of  0.3  [Ref.  1]. 

Several  significant  distinctions  between  static  and 
dynamic  stall  should  be  observed.  Static  stall  is  character¬ 
ized  by  massive  flow  separation  and  loss  of  lift  coinciding 
with  some  critical  angle  of  incidence.  The  aerodynamic 
forces  on  the  airfoil  vary  uniquely,  or  nearly  so,  with  angle 
of  attack.  Dynamic  stall,  securing  during  rapid  increases 
in  incidence,  results  in  a  delay  of  the  associtated  loss  of 
lift  until  well  beyond  the  static  stall  angle.  Dynam.'c  stall 
also  tends  to  be  considerably  more  persistant.  The  flow 
field  does  not  immediately  adjust  to  reductions  in  incidence 
angle  as  with  static  stall.  CRef.  23 

Dynamic  stall  is  characterized  by  the  shedding  and  pas¬ 
sage  of  a  vortex  over  the  upper  surface  of  an  airfoil.  The 
production  of  this  vortex  begins  shortly  after  the  incidence 


has  exceeded  that  of  static  stall.  A  thin  layer  of  reversed 
flow  in  the  boundary  layer,  evolving  into  a  vortex  near  the 
leading  edge,  represents  the  initial  phase  of  a  dynamic  stall 
cycle.  As  this  disturbance  moves  acr:ss  the  upper  surface  of 
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preaches  *:he  trailing  edge,  decrease  dramatically  with  the 
shedding  of  this  vortex.  Reattachment  of  the  flow  will  then 
typically  occur  at  some  angle  less  than  that  corresponding  to 
static  stall.  C  R  e  f  .  23 

B.  THESIS  COALS 

The  research  project  generating  this  thesis  represents  a 
small  cart  of  a  major  investigation  aimed  at  developing  a 
thorough  understanding  of  the  effects  of  compressibility  and 
pitch  oscillation  on  dynamic  stall.  The  primary  goal  of  the 
research  associated  with  thi3  thesis  was  the  development  of 
an  experimental  procedure,  and  supporting  software,  for  the 
acquisition  and  analysis  of  unsteady  flows  with  a  laser  dcp- 
cler  velocimetry  system.  The  procedure  and  software  devel¬ 
oped  was  to  provide  an  initial  analysis  of  possible  areas  of 
interest  and  a  base  upon  which  to  build  as  the  sophistication 
of  the  experimental  methods  increased.  A  successful  attain¬ 
ment  of  this  goal  required  obtaining  velocity  data  in  the 
wake  of  an  oscillating  airfoil  and  wake  profile  plots  pro¬ 
duced  from  this  data. 


I  I  .  EXPERIMENTAL  TONS T  PERAT I  CMS 


A.  EXPERIMENTAL  METHODS 

Since  its  introduction  in  ,Q64,  laser  Doppler  vel- 
:  :  Lnetry  has  proven  to  be  a  highly  reliable  means  for  obtain¬ 
ing  measurements  in  complex  turbulent  flows.  The  nett: i  is 
'apatle  of  performing  instantaneous,  non-in trusive  velocity 
measurements  under  a  wide  range  of  conditions.  Being  rela- 
t  i  v  e  1  y  independent  of  fluid  properties,  accurate  measur  em  er.ts 
are  possible  without-  the  requirement  of  calibration.  IRef.  31 

The  basic  principle  behind  laser  Doppler  velocimetry 
os  the  measurement  of  the  velocity  of  particles  travelling 
with  the  fluid.  These  particles,  when  illuminated  with  a 
focused  laser  beam,  become  sources  of  scattered  light.  The 
Doppler  shift  of  the  scattered  light  is  then  used  to  deter¬ 
mine  the  velocity  of  the  particles  and,  depending  upon  how 


well  *:he  particles  are  following  the  flow,  the  fluid  velocity 


As  a  minimum,  six  components  are  necessary  for  a  com¬ 
plete  LDV  system:  laser  source,  transmitting  optics,  receiv¬ 
ing  optics,  photodetec t or ,  signal  processor,  and  data 
processor.  The  function  of  the  transmitting  and  receiving 
optics  are  focusing  and  light  collecting,  respectively.  The 
photo  detec  tor  serves  to  convert  light  signal  to  an  electrical 


signal  which  is  then  used  by  the  signal  process or  to  provide 
a  measurement  of  the  frequency.  The  capability  of  detecting 
reversed  flows  requires  an  additional  component  for  frequency 
shifting. 

Cue  to  the  wide  range  of  disciplines  involved,  op¬ 
tics,  electronics,  light  scattering,  and  signal  processing, 
*:he  implementation  of  this  method  can  be  difficult  CRef.  33. 
7 onsequent 1 y ,  the  expertise  of  the  operator  can  have  a  pro¬ 
found  impact  on  the  results  of  LDV  measurements .  It  is 
therefcre  imperative.  to  an  even  greater  degree  than  many 
c  the-'  methods,  that  considerable  attention  be  given  to  ap¬ 
plying  correct  measurement  technique. 

2  .  Schl  ierer.  FI  yw  Visualization 

Schliaren  flow  visualization  systems  provide  a  quali¬ 
tative  analysis  of  the  density  gradient  in  compressible  flew 
fields.  'phe  fundamental  concept  of  this  method  is  the  vari¬ 
ation  of  the  refractive  index  gradient  of  a  fluid  with 
density. 

A  basic  schlieren  system  includes  a  light  source, 
transmitting  and  receiving  mirrors,  a  knife  edge,  and  a 
focusing  screen.  The  transmitting  mirror  serves  to  direct  a 
parallel  light  be  am  through  the  test  section  and  onto  the 
receiving  mirror.  An  image  of  the  light  source  is  then 
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therefore,  produce  variations  of  light  intensity  and  provide 
a  visual  representat ion  of  the  flow  field.  [Ref.  41 

B.  INSTRUMENTATION  AND  EQUIPMENT 
1  .  Equipment  Setup 

The  experimental  procedures  for  this  research  were 
conducted  at  the  NASA  Ames  Fluid  Mechanics  Laboratory.  A 
25-cm  by  1 3-crn  indraft  wind  tunnel  with  a  slotted  wall  test 
section  was  used.  Figure  2.1  shows  the  test  section  with 
plenum  chambers  above  and  below  the  slotted  walls.  The 
adaptive  wall  design  of  this  test  section  was  developed  for 
use  in  previous  experiments.  Although  the  technique  was  not 


I  ~ 


t  M  .  *  ' 
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-r 


used,  the  slotted  walls  were  helpful  in  reducing  wall 
interference. 

The  model  for  the  experiment  was  a  NACA  0012  airfril 
j  i  t h  ?”  chord.  The  airfoil,  attached  to  the  test  section 
sails  a*-  quarter  *  h  o  r  d .  was  connected  with  mechanical  linkage 


c  a  i. 

? 1  e  r  *:  r  i  :  n 

-  — 

-  •  » 

as  shown  in  Figure  2 . 1  . 

This  variable 

peed 

m  •:  t  r-  r 
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i  ded 

a  means  for  cscillating  the 

model  about 

P 

i  t  t  h 

axis 

T1 

c: 

obtaining  unsteady  data 

at  varic..s 

r  r  e  q  u  e  r.  c  ies  . 

An  optical  encoder,  also  driven  by  the  motor,  pro¬ 
vided  3,600  counts  per  revolution  of  1  4  -  b  i  t  data  on  the  phase 
jr.gle  of  the  oscillating  airfoil.  The  characteristics  of  the 
:ncoder  were  absolute  in  that  a  whole  word  output,  with  a 
jr.ique  code  pattern,  represented  each  of  the  3,600  discrete 
positions.  This  code  was  derived  from  independent  tracks  on 
:he  encoder  disk  corresponding  to  individual  photodetectors. 

Test  section  Mach  number  and  velocity  were  obtained 
:  y  -using  a  24-port  Scanivalve  with  a  differential  pressure 
rransducer.  Calculation  of  these  values  required  meas- 
.rements  from  four  of  the  ports,  two  for  Mach  number  computa- 
:  i o r.  and  two  for  calibration  of  the  transducer.  The  linear 


-elationship  between  the  measured  pressure  and  the  c  o  .-re¬ 
sponding  voltage  produced  by  the  transducer  was  used  f:r 
let- mining  the  pressures.  In  accomplishing  this,  one  port 
:f  lhe  Coani  valve,  vented  to  the  atmosphere,  measured  ambient 
pressure  and  a  second  port,  “rapped  to  a  mechanical  vacuum 
:umF,  measured  a  calibration  pressure.  A  conversion  factor 


was  then  established  by  dividing  the  differer.  re  of  the  know- 
pressures,  atmospheric  and  calibration,  by  the  difference  in 
voltages  produced  by  the  transducer. 

Location  of  the  LDV  system  probe  volume,  or  beam 
crossing  point.  was  controlled  by  a  three  axis  traversing 
mechanism.  A  commit  Taskmaster  numerical  control  computer, 
originally  developed  for  controlling  movements  of  industrial 
machine  k  cols  ,  was  used  for  this  task.  The  numerical  control 
computer  '  NT?  was  capable  of  manual  operation,  programmed 
operation,  and,  with  software  adaptation,  operation  from  the 
computer  workstation.  The  NCC  provided  control  inputs  to  two 
s?ks  of  electric  motors,  one  for  the  transmitting  optics  and 
the  other  for  the  receiving  optics.  The  necessity  of  using 
different  traversing  motors  for  the  transmitting  and  receiv¬ 
ing  optics  placed  a  great  demand  on  the  synchronisation  cf 
this  system.  Due  to  the  high  sensitivity  of  the  LDV  optics 
train,  the  synchronized  movement  was  often  inadequate  f cr 
maintaining  signal  strength. 

Phase  encoder  data  and  digitized  velocity  data  from 
the  LDV  signal  processors  were  channeled  to  a  3D-LDV  Computer 
Interface  <  C I  5 .  The  Cl,  multiplexing  the  inputs  to  a  single 
digital  data  channel  output,  provided  16-bit  parallel  output 
data  to  the  computer  data  bus.  The  Cl  was  capable  of  accept¬ 
ing  up  to  six  channels  of  digital  data  and  from  one  to  three 
event  pulses  from  the  signal  processors.  The  multiplexing 
cycle  could  be  initiated  by  random  event  pulse  inputs  or 
synchronized  inputs,  the  latter  being  necessary  for  unsteady 


Figure  2.3  LDV  receiving  optics. 


multiline  mode  and  produced  a  vertically  polarized  beam  com¬ 
posed  of  all  wavelengths  inherent  to  this  type  of  laser. 

Figure  2.4  illustrates  the  sequence  of  steps  required 
of  the  transmitting  optics.  A  set  of  mirrors  was  first  used 
to  lower  the  multiline  beam  to  the  level  of  the  optical 
components.  After  passage  through  a  collimator,  used  to 
ensure  collimation  of  the  beam  and  subsequent  crossing  of  the 
beams  at  their  waists,  the  beams  pass  through  an  attenuator. 
The  function  of  this  component  was  to  reduce  the  beam  intens¬ 
ity  without  inhibiting  transmission  of  any  wavelengths,  a 
possible  consequence  of  simply  reducing  the  power  output  of 
the  laser.  Prior  to  separation  of  the  wavelengths  by  the 
dispersion  prism,  it  was  necessary  to  horizontally  polarize 
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the  Bragg  Cell,  the  frequency-shifted  he an  was  deflected  f r : n 
its  initial  alignment.  Therefore,  bean  alignment  prisms,  :r 
steering  wedges,  within  the  Bragg  Cell  nodules  served 
bring  the  be  aims  back  to  parallel. 

A  microscope  : b  :  e  c  t  i  v  e  >  placed  in  the  test  section  at 
"he  probe  volume,  was  used  " o  observe  "he  beans  as  adjust¬ 
ments  were  made  "  o  er.3u.-e  the  beans  focused  to  a  common 
point.  If  "he  two  probe  volumes  are  r.  ot  at  the  same  location 
w  i  "  h  i  r.  "he  tes"  section.  "her.  the  mease,  remerts  will  be  made 
at  do  f ferer"  points.  The  necessity  for  all  four  teams  to 
overlap  at  the  focal  point  was  especially  prevalent  for  the 
acquisition  of  unsteady  data  due  to  the  need  for  a  high 
coincidence  data  rate. 

Receiving  optics  are  designed  to  image,  to  the  great¬ 
est  extent  possible,  the  light  scattered  at  the  probe  volume 
onto  the  photodetectors.  Since  the  scattered  light  intensity 
is  the  greatest  in  the  forward  direction,  a  forward  scatter 
sys"em  is  most  often  preferred  for  obtaining  the  grea"est 
s  i  gr.  a  1  -  to- no  i  se  ratio.  The  system  used  for  this  analysis  was 
an  :ff-axis,  forward  scatter  configuration .  Off-axis  col¬ 
lection  angles  usually  allow  good  noise  isolation  by  a  dm  i  t - 
tir.g  only  light  scattered  at  the  sample  volume  and  improve 
m e a s u r em ent  resolution  by  reducing  the  probe  volume  length. 

;psf  .  i 

The  first  three  lenses  in  the  receiving  optics  were 
used  to  collect  scattered  light  from  the  probe  volume  and 
collimate  the  beam  (Figure  2.5).  The  beam  was  then  direo"ed 


a  polarization  beamsplitter  for 


colors.  The  two  ocmpcnen t s  were  then  passed 
filters,  for  blocking  ambient  light  and  pr- 


talk,  and  to  lenses  used  for  focus!.-.  g 


:  h  e  t  e  am  5 


hole  aperature  of  their  respeoti-re  photomultiplier  tubes. 
The  photomultiplier  tubes  converted  the  light  energy  into 
electronic  signals  which  were  then  passed  to  signal  proces¬ 
sing  equipment.  Preamplifiers  were  located  in  close  proxim¬ 


ity  to  the  photomultiplier  tubes  and  used  to 


mease  the 


signal  level  before  the  introduction  of  additional  noise  into 
the  cables. 

Signal  p recessing  was  accomplished  with  Macrodyne 


ccurter  processors. 


Counter  processors  operate  by  employing 


a  crystal  controlled  oscillator  to  time  over  a  selected  num- 


-f  cycles  within  the  burst  signal.  The  clock  coun* 


initiated  when  the  filtered  signal  exceeds  a  set 


;res  h : . d . 


The  Soppier  frequency  is  then  determined  by  the  time  required 
tc  count  a  specified  number  of  cycles.  The  Macrodyne  proces¬ 
sors  use  a  125  MHz  clock  and  count  over  8  or  1  f>  cycles,  as 


The  Macrodyne  processors  use  a  periodicity  check 
:  u pi e  d  with  a  three-level  validation  scheme  to  reject  noise. 

. e  periodicity  check  uses  two  counters  to  make  period  me  as- 


■ements  over  a  differer. *■  number  of  cycles 


and  8  in  this 


•e.  If  the  raf  i  :  of  the  two  measur  em  ents  is  not  with- 


a  specified  percentage 


-  C  CT  '  O  *-  U 


he  signal  is  rejected 


:  ree- level  validation  requires  a  signal  to  exceed  a  f  c  s  i  t  i  • 


'-\V  v; 

,u  a*  Va  .v*  j 


'-'/v'v'v'v- 

m  I.  *  a >U  fLkJl dr*M  X* JUkJLk. fta ftalUAaiUUUA 


V  *> 


PHOTOMULTIPLIER  POLAR 1 2AT I  ON  COLLIMATING 

tube  beamsplitter  lens 


focusing  lens 

LIME  FILTERS 


FOCUSING  LENS 


HICH 

VOLTAGE 

SOURCE 


PHOTOMULTIPLIER 

TUBE 


HICH  PASS 
FILTER 


w 


PREAMPLIFIERS 


LOW  PASS 
FILTER 


PASSIVE 

V  l 

nil 

JVb. 

(ER 

SINE  WAVE 
GENERATOR 


{  SIGNAL 
|  PROCESSOR 


.[COMPUTER 

•INTERFACE 


COMPUTER 


PROBE 

VOLUME 


G 


OSCILLOSCOPE 


ass 


and  then  exceed  a  neg  a  - 


cld  level,  pass  thrrugh  zero, 
hreshcld  level  before  the  next  zero  crossing.  If  these 
ions  are  not  .net  on  any  cycle  the  counter  will  reset 
.-..red  lately  process  a  new  Bcppler  burst.  This  method 
r  e ; e  c  t  i  o  n  of  pct=-tially  good  signals  having  noisy 
eg  inning  of  the  burst. 


"he  Schliere.n  syst e.m  included  two  concave  mirrors, 
with.  a  '  2  ?  inch  focal  length,  a  stroboscopic  light 
.  a  vertical  knife  edge,  a  focusing  lens,  and  a  polar- 
1m  supcor*-  for  obtaining  photographs  of  the  flow  field, 
rotcsccpe,  a  type  3015  Strobrite,  has  a  maximum  output 
ity  -f  million  beam  candles  and,  using  an  internal 

ator,  car.  provide  flashing  rates  in  excess  of  1,000  per 
.  By  adapting  a  lens  cover  with  a  vertical  slit  to  the 
scope  a  slit  source  was  obtained  and  the  resolution  of 
age  significantly  improved. 

Tc  obtain  the  unsteady  photographs  a  digital  conpar- 
w3s  used  to  trigger  the  stroboscope  and  provide  the 
f:r  the  photographic  imaging  on  the  Polaroid  film.  The 
atsr  received  14-bit  phase  input  from  the  optical  en- 
and  provided  the  triggering  pulse  when  the 
taneous  phase  angle  of  the  airfoil  coincided  with  the 
specified  to  the  c on paratcr. 


-y-V-V  -V -V.v.y. 


DATA  ACQUISITION  AUTOMATION 


A  Digital  VAXstaticn  II  minicomputer  with  graphics  w:rk- 
station  was  used  to  automate  the  data  acquisition  process, 
"he  signal  path  is  illustrated  in  Figure  2.5  starting  from 
the  probe  volume,  through  the  signal  processing  phases,  and 
terminating  within  computer  memory.  The  main  program  used  in 
accomplishing  this  task,  ’LDV_ACQ',  was  originally  developed 
o  o  read  data  from  the  analog  to  digital  converter.  This  was 
later  modified  to  obtain  pressure  data  for  calculation  of  the 
test  section  velocity  and  call  supporting  subroutines  for  the 
acquisition  and  analysis  of  the  velocity  data. 

The  subroutines  assisting  with  various  aspects  of  the 
acquisition  process  are  listed  below  and  included  in  the 
appendix  of  this  report. 

<1  )  HWTRVT  -  allows  control  of  the  traversing  mecha-Ls.m 
from  the  computer  terminal  by  sending  character  data 
to  the  interface. 

(2)  HWWA  -  combined  with  an  internally  called  subrou¬ 
tine,  DAMUX,  retrieves  data  from  multiplexer,  inter¬ 
prets  phase  and  velocity  information,  creates  data 
files,  and  generates  histogram  arrays  and  scaling 
parameters . 

f  3  '  PT_PLOT  -  provides  graphical  display  of  computed 
velocity  values  corresponding  to  airfoil  phase 

angles. 

' A  )  PPLOT  -  similar  to  PT_PLOT  but  provides  the  addit¬ 
ional  capability  of  allowing  the  user  to  vary  the 
plotted  velocity  range. 

<  5  )  HIST  -  plots  a  histogram  of  the  percentage  of  data 
points  which  fall  within  specified  velocity  limits. 
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A.  TEST  SECTION  VELOCITY 

The  main  program  for  the  acquisition  processi  'LDV_ACQ*, 
was  responsible  for  the  calculation  of  test  section  velocity. 
In  combination  with  scanivalve  positioning  subroutines,  the 
program  acquired  pressure  transducer  data  from  7  ports.  As 
previously  discussed,  measurements  from  2  of  the  ports  were 
used  for  calibration.  Pressure  data  from  3  ports  was  dis¬ 
carded,  and  the  other  two  provided  the  static  pressure  at  the 
entrance  to  the  test  section  and  the  total  pressure  measured 
at  the  nozzle  inlet.  Assuming  isentropic,  perfect  gas  condi¬ 
tions,  this  supplied  the  necessary  parameters  for  the  Mach 
number  computation  by  employing  the  following  relationship: 

M  =  (((P0/P)<r*,,/T-l)*(2/(>-l))),/, 

with  a  value  of  1  .4  used  for  7  •  The  test  section  velocity 
was  then  calculated  by: 

V  =  H  x  a 

where  a,  the  sonic  velocity,  was  obtained  as  follows: 

0  =  y/  -yRT 

The  temperature  required  for  this  calculation  was  read  from  a 
wall  thermometer  in  the  test  cell  and  entered  when  prompted 
by  the  program . 

The  velocity  values  were  transferred  to  data  files  and 
later  used  to  normalize  the  LDV  velocity  data. 


B.  SCHLIEREN  PHOTOGRAPHY 


A  qualitative  analysis  of  the  steady  and  unsteady  flew 
fields  was  conducted  at  two  Mach  n urn bers,  0.3  and  0.5,  and 
angles  of  incidence  varying  from  0  to  20  degrees.  In  ac¬ 
complishing  this,  test  section  Mach  number  was  set  by  adjust¬ 
ing  the  wind  tunnel  throat  and  verified  with  Scani  valve 
pressures  and  calculations  from  'LDV_ACQ'.  Final  adjustments 
to  the  Stobosccpe  intensity  and  knife  edge  were  then  made  to 
ensure  optimum  contrast. 

C.  WAKE  VELOCITIES 

With  the  desired  wind  *■  ur.r.e  1  conditions  established  and 
the  probe  volume  in  position  for  the  data  point,  ’HWWA’,  a 
subroutine  called  by  the  main  program,  was  employed  for  raw 
data  acquisition  and  conversion  to  velocity  and  phase  infor¬ 
mation.  An  internally  called  subroutine,  'DAMUX',  acquired  4 
or  6  channels  of  data,  as  selected,  from  the  LDV  Computer 
Interface  (Cl).  The  Cl  was  capable  of  transferring  either  2, 
4  ,  6,  or  8  channels  of  data,  with  the  first  two  of  the 
channels  reserved  for  time  and  status  information  generated 
by  the  Cl.  During  the  acquisition  of  steady  state  velocities, 
a  channels  were  selected  to  provide  the  2  velocity  components 
in  addition  to  time  and  status  data,  which  were  discarded 
since  they  were  not  relevant  for  the  experiment.  For  unsteady 
acquisition,  6  channels  of  data  were  transferred  to  accom¬ 
modate  the  additional  data  channel  required  for  phase  input. 


As  previously  mentioned,  the  signal  processors  determined 


the  Doppler  frequency  by  using  the  time  required  to  count  a 
specified  number  of  cyles,  8  for  this  analysis.  This  raw 


data  was  transferred  to  the  computer  in  16  bits,  with  a  10 
bit  mantissa,  a  4  bit  exponent,  and  2  nonessential  bits  used 
for  controlling  data  flow.  Due  to  the  format  of  the  time 
words,  it  was  necessary  to  invert  and  reverse  the  bits  before 
extracting  the  exponent  and  mantissa  values.  The  signal 

i 

period  was  then  calculated  as  follows: 

ts  =  (  1/32XD9D8.  .  .DO  )  x  2 

where  D9D8...D0  represented  the  binary  mantissa  and  EXP  was 
the  range  selected  with  the  processor,  used  to  divide  the 
number  of  clock  counts  in  the  Doppler  burst. 

The  signal  frequency,  the  inverse  of  the  period,  was  ob¬ 


tained  from: 


f  s  (32  *  10001/C2EXP  x  mantissa  1 


The  signal  frequency  was  composed  of  the  Doppler  frequen¬ 
cy,  the  40  MHz  shift  frequency,  and  the  downmixing  frequency. 


This  allowed  the  Doppler  frequency  to  be  deduced  as  follows: 

^4  ^  I  ^  ,hl,l  ^mli 

The  8  cycles  counted  within  the  Doppler  signal  physically 
represent  8  fringes  crossed  by  a  particle  travelling  with  the 
flow.  Consequently,  the  velocity  component  perpendicular  to 
the  fringe  pattern  may  be  calculated  by  the  product  of  the 


Doppler  frequency  and  the  fringe  spacing.  The  equation  for 
the  latter  is: 


where  A  is  the  wavelength  of  the  bean  and  6  is  the  angle  be¬ 
tween  the  intersecting  beans. 

After  obtaining  a  sufficient  number  of  data  points  at  a 

probe  position  <2,000  for  unsteady  flows),  the  velocity  and 

phase  data,  as  well  as  the  probe  position,  test  section 

velocity,  and  mixing  frequencies,  were  stored  in  files  for 

later  analysis.  The  Taskmaster  traverse  controller  was  then 

« 

used  to  move  the  probe  volume  and  data  collection  continued. 

The  data  acquisition  software  included  two  routines  for 
on-line  performance  checks  of  the  data.  *H1ST'  was  used  to 
provide  a  histogram  of  the  percentage  of  data  points  which 
occurred  within  specified  velocity  limits.  Figure  3.1  shows 
a  sample  plot  provided  by  the  subroutine.  This  data  valida¬ 
tion  method  is  most  useful  for  steady  flow  conditions,  in 
which  case  a  high  percentage  of  data  concentrated  around  the 
estimated  velocity  would  be  a  positive  indicator.  For  un¬ 
steady  analysis  ’PT_PLOT'  provided  a  graph  of  computed  ve¬ 
locity  values  corresponding  to  airfoil  phase  angles.  With 
this  information  available,  the  data  quality  could  be 
easily  assessed  and  a  determination  made  regarding  the  ac¬ 
ceptance  of  these  values. 


VERTICAL  VELOCITY  HISTOGRAM 
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Figure  3.1  Sample  histogram. 
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SCHLIEREN  PHOTOGRAPHS 


A  . 

The  Schlieren  photographs  obtained  from  the  analysis  of 
the  steady  flow  field  were  intended  to  provide  information  on 
the  parameters  corresponding  to  early  flow  separation  and  a 
reference  for  comparison  with  unsteady  flows  at,  the  same 
conditions.  Steady  state  separation  was  first  observed  at  a 
Mach  number  of  0.5  and  a  10  degree  angle  of  incidence  (Fig¬ 
ure  4.1  ).  Based  on  this,  the  parameters  for  the  unsteady 


Figure  4.1  Steady  state  Schlieren  photograph; 
mean  ct  =10  degrees. 


Figure  4.3  Unsteady  flow  Schlieren  photograph; 
f  =  25  Hz,  a  =10  degrees. 

oscillation,  a  logical  result  due  to  the  greater  effect  oi 
time  unsteadiness  at  higher  reduced  frequencies. 


B.  DATA  ANALYSIS  SOFTWARE 

The  main  program ,  ’ORGIZE’,  and  its  supporting  subrou¬ 
tines  were  designed  to  provide  wake  profile  plots  from  thi 
data  files  generated  during  acquisition.  To  accomplish  this 
the  data  file  corresponding  to  each  point  in  the  profile  i; 
input  to  the  software.  A  series  of  routines  first  reorgan 
izes  the  data  files  and  discards  statistically  bad  dat 


po  i  nts  , 


these  outside  of  two  standard  deviations  from  th 


".ear..  Due  to  oscillatory  nature  of  the  forcing  function 

producing  the  unsteady  da*ra,  a  sinusoidal  curve  fit  is  then 
introduced  to  allow  evaluation  of  the  first  ha rm onic  of  the 
velocity  distribution  over  t  he  phase.  A  unique  sinusoid  is 
generated  for  each  jata  point  within  the  profile. 

A  miner  degree  of  uncertainty  is  introduced  into  the  ve¬ 
locities  by  the  data  interpolation  and  curve  fitting  rou¬ 
tines  discussed  above.  However,  unlike  many  other  velooity 
m e a s u r em ent  techniques,  LDV  data  is  obtained  rand om 1 y  due  to 
the  random  arrival  of  the  seed  particles  at  the  prove  volume. 
Without  this  curve  fit,  it  would  be  difficult  to  generate 
velocity  information  at  the  same  desired  phase  angle  for  each 
data  point  in  the  profile. 

During  the  generation  of  the  wake  profile  the  user  is 
provided  a  comparison  plot  of  mean  values  and  the  sinusoid 
curve  fit.  The  mean  values  are  provided  by  ’BDPTS’  and  com¬ 
puted  within  each  20  degree  phase  window,  as  illustrated  in 
Figure  A. 4.  This  was  helpful  in  ascertaining  the  quality  of 
the  curve  fit  and,  consequently,  the  quality  of  the  wake  pro¬ 
file  results . 

The  data  analysis  routines  are  summarized  below  and  in¬ 
cluded  in  the  appendix  of  this  report . 

<  1  )  0RCI2E  -  main  program ,  calls  supporting  subroutines; 
sequentially  organizes  data  to  assist  with  subse¬ 
quent  analysis. 

' 2  '  PT_PL0T  -  provides  graph  of  computed  velocity  values 
c or respond i n g  to  phase. 

'  3  )  BDPTS  -  computes  mean  and  standard  deviation  of  ve¬ 
locities  occurring  within  20  degree  phase  windows! 

^>3 


VLLuUT Y  (ft/s 


HORIZONTAL  VELOCITIES 


discards  data  points  outside  of  two  standard  devia¬ 
tions  of  the  mean  within  each  of  these  windows;  pro¬ 
vides  standard  deviation  plots  of  data. 

f)  A VC  -  averages  data  values  occurring  at  same  dis¬ 
crete  phase. 

S>  INTERP  -  provides  a  linear  interpolation  between 
data »  resulting  in  180  equally  spaced  points;  neces¬ 
sary  for  input  to  FOURIER. 

>>  FOURIER  -  computes  coefficients  for  a  sinusoidal 
curve  fit  of  the  data . (  adapted  from  June  1977  OCS 
Airfoil  program). 

CVFIT  -  fits  sinusoidal  curve  to  velocity  data  ( u- 
tilizes  FOURIER). 


'8)  PLOT  -  provides  graphical  comparison  of  curve  fit 
with  mean  values  generated  by  BDPTS . 

PROPLOT  -  provides  wake  profile  plots. 

WAKE  PROFILES 

L?V  lata  was  obtained  for  steady  and  unsteady  flows  at  a 
Mach  number  of  0.5  and  a  10  degree  angle  of  incidence.  The 
steady  profiles  were  collected  at  distances  of  0.5,  i .0,  and 
2.0  inches  behind  the  trailing  edge  of  the  airfoil.  These 
values  were  non-dimensionalized  by  the  chord  length  of  the 
airfoil  and  are  indicated  as  such  on  the  profile  plots.  Due 
to  time  limitations,  only  one  data  profile  was  collected  for 
the  unsteady  case,  at  two  inches  behind  the  trailing  edge. 

The  wake  profile  plots  are  presented  in  the  following 
pages  of  this  section.  A  lack  of  data  points  in  the  lower 
part  of  the  flow  field  will  be  apparent  upon  inspection  of 
these  results.  This  was  due  to  a  mechanical  restriction  on 
*rhe  traversing  mechanism. 

The  curve  fit  of  the  plotted  values  of  the  steady  data 
was  approximated  and  hand  drawn  to  aid  the  reader  in  ob¬ 
serving  the  effects  of  the  airfoil  in  the  flow  field.  The 
software  provides  a  spline  interpolation  fit  between  the 
data,  as  used  in  the  unsteady  plots.  However,  due  to  the 
amount  of  data  scatter  involved,  this  method  was  not  used  for 


these  results. 
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Figure  4.8  Steady  wake  profile. 
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Figure  4.12  Unsteady  wake  profile, 
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The  purpose  of  the  results  obtained  during  the  experi¬ 
mental  phase  of  this  thesis  was  to  validate  the  technique  and 
software  developed.  This  was  accomplished  when  unsteady 
data,  generated  by  an  oscillating  airfoil,  was  successfully 
acquired,  analyzed,  and  reduced  to  wake  profile  plots. 

Throughout  the  course  of  this  investigation  hardware 
problems  have  a  posed  a  signigicant  challenge  to  acquiring 
data  within  the  existing  time  constraints.  Consequently,  the 
data  obtained  was  accepted  and  included  in  this  report  with¬ 
out  the  opportunity  to  more  carefully  investigate  regions  of 
significant  data  scatter  and  improve,  or  validate,  the  data 
plots  . 

'rhe  results  produced  by  these  procedures  confirm  a  sound 
technique,  but  also  reveal  several  areas  in  need  of  addi¬ 
tional  emphasis.  At  present,  an  interpolation  subroutine  is 
required  to  produce  data  at  equally  spaced  intervals  before 
generation  of  the  curve  fitting  coefficients  by  ’FOURIER'. 
Combined  with  the  inability  to  utilize  a  large  number  of  data 
points  (  <  300),  this  limitation  will  inherently  introduce 
additional  uncertainty  into  the  results.  A  second  software 
modification  to  be  considered  is  the  generation  of  the  veloc¬ 
ity  data  at  various  phase  angles  simultaneously.  Currently, 
the  data  must  be  processed  for  each  separate  phase  angle 


desired.  This  would  eliminate  the  need  for  storage  of  large 
amounts  of  raw  data  and  significantly  reduce  the  processing 
time.  By  increasing  the  range  of  vertical  mov em ent  of  the 
traverse,  the  mechanical  restriction  that  prevented  the  ac¬ 
quisition  of  data  in  the  lower  regions  of  the  wake  can  be 
el  iminated  . 

Continuation  of  this  study  will  be  accomplished  with  a 
newly  developed  test  section  currently  being  installed.  This 
test  section  will  enable  schlieren  flow  visualization  and  LDV 
data  acquisition  at  locations  over  the  airfoil  which  were 
previously  blocked  by  the  mechanical  linkage  used  for  sup¬ 
porting  and  oscillating  the  airfoil.  This,  combined  with  the 
modifications  suggested  above,  will  enable  an  in-depth  analy¬ 
sis  of  wake  and  boundary  layer  characteristics  associated 
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APPENDIX 


COMPUTER  PROGRAMS 


PROGRAM  LDV  ACO 


Adapted  from  LAB STAR  Exaaple3_l  by  S.  D.  Hedrick  in  Hay,  1987  (or  reading 
A-D  data.  Hodified  in  AuguatT  1987  to  be  uaed  at  the  central  prograa  in 
the  acquisition  of  uneteady  velocity  coaponenta  with  an  LDV. 


INCLUDE  'eye$library:LIOSET. FOR’ 

C  declare  local  variablea 

INTEGER  sysstat  (status  ret  by  LIO  calla 

INTEGER  dev_id,  elk  id  ILIO  device  ID 
INTEGER  data_length“  inumber  of  bytea  of  data  read 
INTEGER  i,  npt,  dnpt 
INTEGER  ehn,  natepa,  ier 

REAL  rate,  volt(0:23>,  press(0:23),  corfac 

REAL  aua,  aean ,  calp,  conv,  aab,  teap,  aach,  velo 

CHARACTER* 1  ana 

CHARACTER  *  4  clock 

LOGICAL  hoae 


C  declare  data  buffers 

I NTEGER • 2  raw  data(SOl) 
REAL  VOLTAGES (501) 


OPEN! UNIT-1  3 , riLE-' scan i  .dat'  , STATUS- ' NEW’  ) 

write!*,*)  'ENTER  A-D  CHANNEL  TO  BE  USED  (0,1,8.91' 
read!  • ,  * )  chn 

10  write!*,*)  'ENTER  CLOCK  RATE  (.01  -  4500  Hz ) --ALLOW  TOR  A', 

1  '  HAXIHUH  Or  2  SECONDS  TO  ACQUIRE  DESIRED  SAMPLES' 

read! • , •  )  rate 

if  (rate  .It.  .01  .or.  rate  .gt.  4500)  then 

write!*, *1  'RATE  IS  NOT  WITHIN  SPECIFIED  LIHITS' 
go  to  10 
end  i  f 

20  write! *,*) 'ENTER  NUMBER  or  DATA  POINTS  TO  BE  COLLECTED', 

1  '  (  500  aax )  ’ 

read! •  ,  •  )  npt 

if  (npt  .It.  1  .or.  npt  .gt.  500)  then 

write!*,*)  ' NUHBER  OR  POINTS  IS  NOT  WITHIN  LIMITS' 
go  to  20 
end  i  f 

npt  -  npt  ♦  1 
dnpt  -  npt  •  2 

write!*,*)  'ENTER  CLOCK  DEVICE  (KZAO.KZBO)’ 
read! • , ' ( A4 ) ' )  clock 


Attach  the  ADVU-D  and  set  up  for  queued  I/O 

ays  a ta t-L IOSATTACH (dev  td.'AZAO',  LIOSK  qio) 

I r ( .NOT.  sys_stat)  CALL  LIBSSIGNAL! kval ( sys_stat ) I 


Set  up  the  ADV  for  synchronous  transfer  on  selected  channel 
Gain  of  1 


ays  atat  -  LIOSSET  I  (dev  id,  LIOSK  SYNCH ,  0) 

If! .NOT.  sys_statT  CALL  LIBSSIGNAL! kval  (  ayaatat ) ) 


46 


non  non  noon  non  onooo 


sys_stat  •  L I  OS  S£T_ I  (devoid,  LIOSF  AD_CHAN ,  1,  chn) 

in. NOT.  sys_stat7  CALL  LIBSSIGNALI tval  i  sys_stat )  ) 

SyS  Stat  -  LIOSSET_I  (dev  id,  LIOSF  ADGAIN,  1,  1) 

in. NOT.  sys_statT  CALL  L I B  JS I GNAL  (  l  va  1  (  sys_s  ta  t )  ) 

Get  a  raw_data  buffet 

length  of  buffer  and  returned  data  length  is  in  bytes 
Attach  clock  and  set  rate 

home  -  .true, 
nsteps  -  0 
do  40  }-0 , 7 

call  dal  1 j 1 lab( home , ns teps , ie r ) 

sys  stat-  L I  OS ATTACH (clk_id, clock, LI OSK_QI 0 ) 

in  -NOT.  sys  stat)  CALL  LIBSSIGNALI tval ( sysstat ) ) 
sys  stat-LIOSSET  R( elk_id, LIOSR_CLK  RATE,  l,rate) 

in .NOT.  sys  stat)  call  LIBSSIGNALI Ival I sys_stat ) ) 
sys  stat-LIOSSET  Ilclk  id , LIOSR_TRIG ,  1,  LIOSK  IMMEDIATE ) 
in. NOT.  sys  stat)  CALL  LIBSSIGNALI  %vaT(  sys_stat )  ) 
sys  stat  •  L I  OS  READ ( 3e v  id.  raw  data , dnpt , da ta  length,  ) 

in.  NOT  .  sys_stat)  CALL'LIBSSIGNALI  kvaT  (  sys  stat)) 
sys  stat  -  LIOSSET  I  (dev_id,  LIOSK_TRIG,  2,  LIOSR_£lK_POINT, 
1  elk  id! 

I r (  .NOT?  sys  stat )  CALL  L I BS S I GNAL ( \ va 1 ( sy S_S t a t )  ) 
Convert  the  raw  data  to  voltages 

call  LSPSFORMAT  TRANSLATE_ADC( raw_data, voltages, npt,  ,  ,  ) 


calculate  the  mean  value  of  the  voltages  and  detach  the  clock 

sum  -  0.0 
do  30  i-2, npt 

sum  -  sum  ♦  voltagesli) 

30  continue 

mean  -  sum  /  float(npt-l) 
vol t ( j )  •  mean 

sysstat  •  LIOSDETACHI elk  id,  ) 

IF!  .NOT.  sys_s tatT  CALL  L I BS S I GNAL I t va 1 1 sy s_s t a t > ) 

40  continue 

detach  from  the  A/D 

sys  stat  -  LIOSDETACHI dev  id,  ) 

IE! .NOT.  sys  statT  CALL  LIBSSIGNALI %val ( sysstat ) ) 

convert  scanivalve  voltages  into  units  of  pressure 

write!*,*)  'ENTER  CALIBRATION  PRESSURE  (in.  of  Hg)' 
read! * , * )  calp 

write!*,*)  'ENTER  REQUIRED  CORRECTION  TACTOR  FOR  CALIBRATION' 
1  '  PRESSURE  (in.  of  Hg)' 

read!*,*)  corfae 
calp-calp*corfac 

write!*,*)  'ENTER  ATMOSPHERIC  PRESSURE  (in.  of  Hg)' 
read!*,*)  a mb 
press(O)  -  calp 


non 


press ( 1 )  -  sab 

conv  -  t  press ( 1 ) -press !  0 ) )  '  ( vo 1 1 ( 1 ) - vol t ( 0 ) ) 

do  50  , 7 

press!])  *  conv  •  volttj)  «  calp  -  voltlO)  •  conv 
50  continue 

WRITE!  •  ,  *  I  VOLT!  0 )  ,VOLT( 1 )  , VOLT! 2 )  , VOLT! 7 ) 
wr  i  te( 1 3 , 55 ) 

5  5  format!  2X,  '  PORT'  , 4X, ' PRESSURE' ,//) 

do  60  i«0  ,  7 

write! 13.  '  < 3X. 12, SX, T9 . 4) ’ )  i,  press(i) 

60  continue 

coapute  test  section  velocity 

write! *,•) 'ENTER  AMBIENT  TEMPERATURE  (deg.  T ) ' 
read! • , * )  teap 

aach«sqrt(abs( 5.*! (press! 2)/press!7) )••(  .  4/1 . 4  )  -1 .  )  )  ) 
vel o»aach*sqrt(l. 4*53. 3*32. 174*(teap«-459. 67)) 

write! *,») 'TEST  SECTION  MACH  NUMBER  AND  VELOCI TY-- ' , aach , ve 1 O 

wr i te( *,*) 'CONTINUE  WITH  LDV  DATA  ACQUISITION?  ( Y/N ,  DEr-Y)' 
read! • , ' ( A1 ) ' )  ans 

if  (ans  .eq.  'n'  .or.  ans  .eq.  'N')  go  to  70 
C 

C  call  subroutine  to  eove  LDV  traversing  aechanisa  and  acquire  data 

c 


c 


call  HWTRVT(velo) 


SUBROUTINE  HWWA ( XROS , Z  POS , VE  LO ) 


C  acronym :  Hardware  checkout  -  test  of  DFVII-wa 
C  PURPOSE: 

C  r or  LVIS,  hwdvyb  drives  the  acquisition  of  the  LDV  data. 

C  It  is  a  test  program  for  the  YB  driver. 

C 

C  METHOD: 

C 

C  BEGIN 

C  Install  loadable  drivers  if  necessary. 

C  Read  the  various  instruaent  types  in  succession,  skipping  those  which 
C  are  not  to  be  read. 

C  Reaove  loadable  drivers  if  necessary. 

C  END 

C 

C  PARAMETERS: 

C  ARG  DIM  TYPE  I/O/S  DESCRIPTION 

r 

C  COMMONS  USED: 

C 

C  FILES  USED:  None 

C 

C  ERROR  HANDLING: 

C  IER  -  1  aeans  no  "fatal*  error. 

C 

C  NOTES: 


C  LOCAL  VARIABLES: 


r 

r 

VAR 

DIM 

TYPE  DESCRIPTION 

C 

EXTERNAL  REFERENCES: 

r 

NAME 

DESCRIPTION  AND  SOURCE 

n 

r 

DVWA 

Calls  DRV11-WA  device  driver  to  acquire  A/D  data. 

C 

c 

STANDARDS  VIOLATIONS: 

None . 

r' 

r * 

ENVIRONMENT: 

DEC  VAX  /VMS  AND  FORTRAN  77 

DEVELOPMENT  HISTORY: 

r 

DATE 

INITIALS 

DESCRIPTION 

r 

07/26  '83 

THL 

Adapted  froa  DAPNT 

r 

12/12/85 

CLH 

Modified  to  print  results  in  either 

r 

octal  or  integer. 

c 

02/25/87 

GBG 

Modified  for  use  with  DRV11-WA  and  VMS 

c 

05/06/87 

CLH 

Modified  to  accept  "new*  DAMUX1 . 

c 

06/08/87 

RRR 

Modified  to  subroutine  ♦  deleted  CLH  5/87  Mod 

r 

JUL-AUG, 1987 

SDH 

Modified  to  calculate  velocities  using  variable 

mixing  frequency  as  user  input,  interpret  phase 
angle  inforaation  generated  by  digital  encoder, 
transfer  test  section  velocity  and  probe  posit¬ 
ion  to  data  file,  and  generate  histograa  arrays 
and  scaling  paraaeters 


AUTHOR ( S ) :  Ted  Lichtenstein,  inforaatics  General  Corp. 


IMPLICIT  NON* 


> 


•  « 


c 


c 


c 

c 

100 


PARAMETER 

PARAMETER 

PARAMETER 


NRAWMX-80000 

LUNTI-5 

LUNOUT-6 


I NTEGER  *  2 
I NTEGER • 2 
I NTEGER  *  2 
INTEGER 
INTEGER 
INTEGER 
INTEGER 
INTEGER 
INTEGER 
REAL 
LOGICAL 
CHARACTER  *10 
CHARACTER  *  1 
REAL*  4 
REAL*  4 
INTEGER 
REAL 


I RAWDT! NRAWMX) ,  I  FORK ,  MUXTIME  (16  384) 

MUXSTS  ( 16384 ) , RAWDAl 20000 ) , RAWDBI 20000 ) 

RAWDC ( 20000 ) 

STATUS, I ER , PHAN 

NCOLUM,  IYESNO.  NRAWIN,I,J 

LIBSSUBX,  NCHANS,  NSAMPLES ,  INDEX,  JNDEX 

MRAWA , MRANB , MRAWC , MRAWD 

I RAWDTA , I RAWDTB , I RAWDTC , I RAWDTD 

IRAWDTE, IRAWDTT 

PREQ( 2) , VEL ( 2 ) , MI XA ,KIXB,XPOS, ZPOS , VELO 
DATSELCT  ( 3 ) 

INPILE 

ANS 

yu( 160) , rv( 120) , yui ( 160 ) ,rvi ( 120 ) , ymaxu, ymaxv 

SAVE 

COUNT , XLOW , XH I 
X(  4000 ) , Y ( 4  0  0  0 ) 


INCLUDE  ’ ( SSYSSRVNAM ) ' 

WRITE! *,*)’ ENTER  DESIRED  NAME  OP  DATA  PILE--E.G. ,  LDVXXX.DAT,’, 

1  '  WHERE  THE  ‘XXX‘  MAY  BE  USED  AS  A  SERIAL  ID  POR  THE', 

1  '  PROBE  POSITION' 

READ! • , ' ( A10 ) ' )  INFILE 

open  (unit-11, file-INPILE,status-'new') 


CONTINUE 

C  •  Fill  the  raw  data  buffer  with  this  bit  pattern: 
C  •  ‘1010101010101010* 

C  *  This  equals  125252(8),  -21846(10). 

C 

DO  200  I  -  1,  NRAWMX 

IRAWDTII)  -  ' 1 25252 ’ O 
200  CONTINUE 


DO  I  -  1 ,  16384 

MUXTIME  (I)  -  '  1  252  52  ' O 
MUXSTS  (I)  -  '  125252 ' O 


END  DO 


C 

220 

C  •  Proapt 

C 


1 

1 


c 


CONTINUE 

for  no.  of  words  to  acquire: 

WRITE!*,*)'  ' 

WRITE! *,*)' Starting  LDV  acquisition' 

WRITE!  *,*)'  ' 

WRITE! *,*)' ENTER  NUMBER  OP  CHANNELS  TO  BE  USED', 

'  (assumes  non-osci 1 la t i ng  airfoil  if', 
'  lees  than  six)' 

READ!*,*)  NCHANS 


50 


Fvv 


WRITE  (*,  225) 

FORMAT  ('  Enter  number  of  samples:') 

READ  I*.  •)  NSAMPLES 

ir  (NSAMPLES  .GT.  20000)  THEN 
WRITE  (*,  226) 

FORMAT  I'  The  number  of  samples  must  be  20,000  or  less.' 
GO  TO  224 

END  ir 


CONTINUE 

Prompt  for  output  format  of  data: 

I  FORM-2 

IT  (  IFORM.LT. 1  .OR.  ITORM.GT.2  )  GO  TO  230 


Select  the  data  to  be  returned. 


DATSELCT  (1)  -  .false. 

DATSELCT  (2)  -  .false. 

DATSELCT  (3)  -  .TRUE. 

Prompt  for  the  mixing  frequencies  selected  by  user 

WRITE) «,♦>' ENTER  MIXING  FREQUENCIES  FOR  EACH  COMPONENT  -  (blue, 
lgreen)  -  in  MHz' 

READ ( * , * )  HIXA.MIXB 


Read  and  display  the  data  one  sample  at  a  time. 

IF  (NCHANS  . LT.  6)  THEN 

WRITEt *,*)’ ENTER  AIRFOIL  PHASE  ANGLE  (tenths  of  degrees,  14)’ 

READ (  •  ,  *  I  PHAN 
DO  9  I  - 1 , NSAMPLES 

CALL  DAMUX (NCHANS, 1 .DATSELCT, MUXTIME, MUXSTS, I RAWDT, I ER) 

RAWDA (  I  ) -  I RAWDT  (  1 ) 

RAWDB (  I  )  - I RAWDT  (  2 i 
RAWDC ( I ) - PHAN 
CONTINUE 
GO  TO  15 
END  IF 

DO  10  1-1, NSAMPLES 

CALL  DAMUX  (NCHANS,  1.  DATSELCT,  MUXTIME,  MUXSTS, 
•  I RAWDT,  I ER ) 

RAWDA ( I  1  —  I RAWDT (  1  ) 

RAWDB ( I ) - I RAWDT ( 2 ) 

RAWDC ( I ) -I RAWDT! 3 ) 

CONTINUE 

DO  19  1-1, NSAMPLES 
MRAWA-0 


19 

15 


HRAWB-0 
HRAWC-0 
HRAWD-0 

HRAWA- 1 1 B I TS ( RAWDC (  I  )  , 1 2  ,  2 ) 

HRAWB-I I  BITS! RAWDC 1  1  )  , 8 , 4 ) 

HRAWC-IIBITS! RAWDC ( I ) , 4 , 4 ) 

HRAWD-I I  BITS! RAWDC!  I  )  . 0 , 4 ) 

RAWDC  (  I  )-HRAWA*1000*HRAWB*100*HRAWC*10*HRAWD 
CONTINUE 

DO  20  I  •  l , NSAHPLES 

RAWDA (  I  ) • I NOT ( RAWDA(  I )  ) 

RAWDB  I  I ) • I NOT ( RAWDB I  I )  ) 

IRAWDTA-0 
IRAWDTB-0 
IRAWDTC-0 
IRAWDTD-0 
IRAWDTE-0 
IRAWDTF-0 

I RAWDTA- 1 1 B ITS ( RAWDA (  I  > 

I RAWDTB- 1 1 B I TS ( RAWDA (  I  ) 

I RAWDTC- 2**1 RAWDTA  *  I RAWDTB 
IF  (I RAWDTC  .£0-  0)  THEN 
WRITE! » , * ) ' IRAWDTC-0' 

GO  TO  20 
END  IF 

FREQ ( l)-((32.*10.**3.)/( FLOAT! I RAWDTC ) ) )-40 . ♦NIXA 
I RAWDTD- 1 1 B I TS ! RAWDB  (  I  ) , 1 0 , 4  ) 

IRAWDTE-I I  BITS! RAWDB!  I  )  , 0 , 10 ) 
IRAWDTF-2**IRAWDTD*IRAWDTE 
IF  (IRAWDTF  .EQ.  01  THEN 
WRITE! *,«) ' I RAWDTF-0 ' 

GO  TO  20 
END  IF 

FREQ! 2)-( ( 32.*10.**3. )/( FLOAT! IRAWDTF) >  )-«0 . *HIXB 
vel(l)-freq(2)*0. 5145/(2. *sin(atan(13./(2. *482. 6)))) 
vel(2)-£req(l)*0. 4880/(2. *sin(atan(13./!2. *482. 6)))) 
vel ( 1 )«vel ( 1 )/. 3048 
vel(2)-ve 1(2)/. 3048 


,10,4) 

,0,10) 


20 


WRITE! 11 , 
WRITE! 11 , 
WRITE! 11 , 


WRITE! 11 , 700)  vel ( 2 ) , vel ( 1 ) , RAWDC! I ) 
CONTINUE 
900)  XPOS.ZPOS 
950)  VELO , NSAHPLES 
975)  HI XA , N I XB 


WRITE! *,*)' DO  YOU  DESIRE  A  HISTOCRAH  OF  THE  VELOCITIES  OBTAINED? 1 
'  ( Y/N  -  DEF-N) ' 

A1 ) ' )  ANS 

NE.  'y'  .AND.  ANS  . NE .  'Y')  THEN 

5 


40 


READ!  •  ,  '  ( 
IF  (ANS  . 

GO  TO  3 
END  IF 
CLOSE ( UN I 
OPEN! UNIT 
DO  40  I-l 

READ! 11 
CONTINUE 
XLOW--200 
XHI--19S 
DO  80  J-l 


T-ll  ) 

-11 .FILE- INFILE, STATUS- ’OLD’ 

.NSAHPLES 

, 800 )  X! I ) , Y< I ) 


160 


52 


■.  s.-\  f .  . 

.  F  >  •  a.  -  W  r  -  ft.  »  >  '  , 


V?* 


••  .•  .  • » V/’ 

a\»  F-  — .V  mF.  m\  m  A  — %  -V  -"a  AiV  Al>  A*H  alV.M,V 


COUNT*  0 

DO  50  1-1, NSAHPLES 

I r  ( X ( I )  .GE.  XLOW  .AND.  X ( I )  .LT.  XHI)  THEN 
COUN T- COUNT *1 
END  I  F 
50  CONTINUE 

YU  I J  i  - i  FLOAT! COUNT )  /NSAHPLES ) *100 . 

XLOW-XLOW*  5 
XH I -XH 1+5 
60  CONTINUE 

Calculating  a  acaling  parameter  for  the  horizontal  velocity  histogram 

DO  62  1-1.160 
YU  1 (  I  ) -  YU (  I ) 

62  CONTINUE 

DO  63  1-1,160 
DO  65  J-l  .  159 

IT  (YUKJ  +  l)  .LT.  YU1(J>>  THEN 
SAVE- YU1  (  J  ) 

YU1 l J I -YU1 ( J*1 ) 

YU1 ( J*1 I -SAVE 
END  IT 
65  CONTINUE 

63  CONTINUE 
YHAXU-YU1  ( 160 1*10.0 

XLOW--300 
XHI--295 
DO  80  J-l, 120 
COUNT-0 

DO  70  1-1, NSAHPLES 

IT  (Y(I>  .GE.  XLOW  .AND.  YII)  . LT .  XHI)  THEN 
COUNT-COUNT* 1 
END  ir 
70  CONTINUE 

YV< J  >•< FLOAT ( COUNT ) /NSAHPLES 1*100. 

XLOW-XLOW*  5 
XHI -XH 1*5 
80  CONTINUE 

Calculating  a  scaling  parameter  for  the  vertical  velocity  histogram 

DO  82  1-1,120 
YV1 ( I ) - YV ( I ) 

82  CONTINUE 

DO  83  1-1,120 
DO  85  J-l  ,  119 

IT  (  YV1  (  J*1 )  .LT.  YVKJ!)  THEN 
SAVE- YV1 ( J ) 

YV1 ( J )-YVl ( J*1 ) 

YV1 ( J*1 ) -SAVE 
END  IF 
85  CONTINUE 

83  CONTINUE 
YHAJCV-YV1  (  120  1*10 . 0 

CALL  HI ST( YU , YV , YHAXU , YHAXV ) 


35  CLOSE! UNIT-1 1 ) 


WRITE! *.*)' DO  YOU  DESIRE  A  PLOT  OP  THE  COMPUTED  VELOCITY', 

1  '  VALUES?  ( Y/N  -  DEr-Y)’ 

READ! • , ' ( A1 ) ’ I  ANS 

IP  ! ANS  . EQ .  'N'  .OR.  ANS  . EQ .  'n')  THEN 
GO  TO  30 
END  I  P 

CALL  PT_PLOT( NSANPLES , INPILE) 

WRITE! • , • ) ' DO  YOU  DESIRE  TO  VARY  THE  PLOTTED  VELOCITY  RANGE?', 
I  '  ( Y/N--DEF-Y ) ' 

READ' • , ’ ( A1 ) ' )  ANS 

IP  (ANS  . EQ .  'N'  .OR.  ANS  . EQ .  'n')  GO  TO  30 
CALL  PPLOT! NSAMPLES , INPILE ) 

•  Again? 

30  WRITE! • , • ) ' 

WRITE  (LUNTI,*)  'Again?  (O-NO,  1-YES)  ' 

READ  (LUNTI,*)  IYESNO 

IP  (  IYESNO. EQ.l  )  GO  TO  100 

RETURN 

*  For«at* - 

■»0  0  PORNA.  .  3X,  .  9 . 4 , 1 X , F9 . 4 , 8X , I  4 ) 

800  FORMA.  '.  ,X,  F9 . 4  ,  IX  ,  F9  .  4  ) 

900  FORMAT! //, 2X,  ' XPOS- ’ , F6 . 2 , 2X, ' ZPOS- ' , P6 . 2 ) 

950  FORMAT! 2X, ' TEST  SECTION  VELOCITY- ', T5 . 1 , 2X , 15 , '  SAMPLES') 

9’S  FORMAT! 2X,  ' BLUE  AND  GREEN  MIXING  FREQUENC I ES-- '  , 2  PS . 1 ,  ' ( MHz )  ’  ) 
END 


SUBROUTINE  PT  PLOT! NSAMPLES , INFI LE ) 


written  by  S.  D.  Hedrick  in  July,  1987  to  plot  computed  velocity  values 
corresponding  to  airfoil  phase  angles.  This  subroutine  is  intended  to 
provide  the  user  with  graphical  information  on  acquired  LDV  velocities 
and  assist  in  making  a  determination  regarding  the  adequacy  of  the  data 
obtained  at  the  current  probe  position. 


IMPLICIT  NONE 
INTEGER  NSAMPLES,  I 

REAL • 4  XVALl 4  000 )  , YVALU!  4  000  )  , YVALV!  4  000  )  ,XCONTR( 4 )  , YCONTRl 4 ) 

CHARACTER  *  1 0  INPILE 
CHARACTER • 1 1  X LABEL 
CHARACTER* 15  Y LABEL 
CHARACTER *19  TITLE 
CHARACTER* 1  DUMMY 

OPEN  ( UN  IT- 11 , PILE- INTI LE , STATUS- 'OLD'  ) 

DO  10  1-1 , NSAMPLES 

READ  I  11,100  1  YVALU (  I  )  , YVALV ( I  )  , XVAL I  I ) 

XVALl  I  I -XVALl  I  1/10.0 
10  CONTINUE 

)0  PORMATI SX,P9.4,1X,P9.4,8X,P4.0) 

Set  up  the  coordinate  axis 

XLABEL- '  PHASE  ANGLE' 

YLABEL- ' VELOCITY  (ft/S)' 

TITLE- ' COMPUTED  VELOCITIES' 

XCONTRI 11-6.0 
XCONTRI  2  1-0.0 
XCONTRI 31-360.0 
XCONTRI  4  1-45.0 
YCONTRl  11-6.0 
YCONTRl 2 I--200 . 0 
YCONTRl 3 1-600.0 
YCONTRl  4  1-100.0 

CALL  LGPSPLOTI 1 ,' EXSY' , XVAL, YVALU, 0, XLABEL, YLABEL,  ,  ,  ,XCONTR, 

1  YCONTR ,  .TITLE) 


Plot  the  individual  data  points  on  the  axis 

CALL  LGPSPOINTI 1 , XVAL, YVALU, . 03 , 3 , NSAMPLES ,  , 
CALL  LGPSPOINTI 1 , XVAL, YVALV, . 03 , 4 , NSAMPLES ,  , 

CALL  LGP$PUT_TEXT( 1,4.5,5.75,'*  U-COMPONENT' , 
CALL  LGP$PUT_TEXT( 1 , 4 . 5, S. 5, ' 0  V-COMPONENT' , 

WRITE!*,*)  'Press  <CR>  to  terminate  plot' 

READ! *,200)  DUMMY 
200  PORMATI All 

CALL  LGPSTERMINATE  PLOT! 1,1) 

CLOSE! UN IT- 11 ) 

RETURN 

END 


.*  V.  v-  »*  . 


rlvl'-- -IvCvIvIjN. 


o  o  o  non 


SUBROUTINE  PPLOT ( NSAHPLES ,  I  NT  I LE ) 

IMPLICIT  NONE 
INTEGER  NSAHPLES ,  I 
REAL  VM I N  ,  VMAJC  ,  I NTV 

REAL *4  XVALI  4  000 ) . YVALU!  4  0  00  ),  YVALVI 4000 }  , XCONTRI  4 )  , YCONTRI 4 ) 

CHARACTER  *10  INPILE 

CHARACTER  *11  X LABEL 

CHARACTER • 1 S  Y LABEL 

CHARACTER  *19  TITLE 

CHARACTER* 1  DUMMY 

OPEN  ( UNIT-11 , TILE-INFILE , STATUS- ' OLD' 1 
DO  10  1-1, NSAHPLES 

READ (  11,100)  YVALU!  I  ) , YVALV( I )  , XVALI  I  ) 

XVALI  I  l-XVALI  I  1/10.0 
10  CONTINUE 

100  FORMAT! 5X,F9.4,lX,P9.4,8X,r4.0) 

Set  up  the  coordinate  axis 

WRITE! *,*)' ENTER  VELOCITY  RANGE  DESIRED  |lin,aax) ’ 

READ!*,*)  VRIN, VMAX 
INTV- I VMAX-VMIN ) /S . 

XLABEL- ' PHASE  ANGLE' 

YLA8EL- ' VELOC I TY  (ft/S)' 

TITLE- ’ COMPUTED  VELOCITIES' 

XCONTRI  11-6.0 
XCONTRI 21-0.0 
XCONTRI 3 1-360 .0 
XCONTRI 4 1-45.0 
YCONTRI 11-6.0 
YCONTRI 2  I -VM I N 
YCONTRI 3 1 -VMAX 
YCONTRI  4 1  -  I NTV 

CALL  LGPS  PLOT  I  1 ,  ' EXS Y '  , XVAL , YVALU , 0 , XLABEL , YLABEL ,  ,  ,  , XCONTR , 

1  YCONTR ,  .TITLE) 

Plot  the  individual  data  points  on  the  axis 

CALL  LGPSPOINTI 1 ,XVAL, YVALU, .03 , 3 .NSAHPLES,  ,  ,  ) 

CALL  LGPSPOINTI 1 , XVAL, YVALV, . 03, 4 .NSAMPLES,  ,  ,  ) 

CALL  LGPSPUT_TEXT( 1 , . 05, . 3, 'ASTERISK-U  COMPONENT’,  ,  ) 

CALL  LGPSPUT_TEXT( 1 , .05, . 1 , 'CIRCLE-V  COMPONENT’,  ,  ) 

WRITE!*,*)  'Press  <CR>  to  terminate  plot’ 

READ! *,200)  DUMMY 
200  FORMAT (All 

CALL  LGPSTERMI NATE_PLOT  (  1,1) 

CLOSE  I  UN I T-l 1 ) 

RETURN 

END 


non  non 


SUBROUTINE  H I  ST  (  YU ,  YV ,  YMAXU  .  YMAXV  ) 


Written  by  S.  D.  Hedrick  in  August,  1987  to  plot  a  histogram  of  the  per¬ 
centage  of  data  points  which  fall  within  specified  velocity  limits.  The 
subroutine  produces  histograms  for  each  of  two  velocity  components.  Scal¬ 
ing  parameters  for  the  vertical  axis  and  data  arrays  are  computed  within 
a  preceding  subroutine  (HWWA). 


IMPLICIT  NONE 

REAL* 4  YU(160),YV(120) ,XCONTRU( 4 ) ,XCONTRV( 4 ) , YCONTRU! 4 ) 

REAL » 4  YMAXU.YMAXV, YCONTRV!  4  I 

REAL  XLOW , XH I , XLOWU ( 160),XHIU(160) ,XLOWV( 120). XHX V( 120 ) 

INTEGER  I 

CHARACTER  *  29  UTITLE 
CHARACTER* 2 7  VTITLE 
CHARACTER  *10  X LABEL 
CHARACTER*18  YLABEL 
CHARACTER* 1  DUMMY 

Setting  up  the  coordinate  axis 

XLABEL-' VELOCITIES' 

YLABEL-' PERCENTAGE  0 T  DATA' 

UTITLE- ' HORIZONTAL  VELOCITY  HISTOGRAM' 

VTITLE-’ VERTICAL  VELOCITY  HISTOGRAM' 

XCONTRU ( 1 )-6 . 0 
XCONTRU ( 2  I --2  00 . 0 
XCONTRU ( 31-600 . 0 
XCONTRU! 41-100.0 
XCONTRV ( 1 ) -6 . 0 
XCONTRV (  2  )  —  300 . 0 
XCONTRV!  3  ) -  300 . 0 
XCONTRV! 41-100.0 
YCONTRU ( 1 )-6 .0 
YCONTRU  C  2 ) —  0 . 0 
YCONTRU! 3 I -YMAXU 
YCONTRU! 41-10.0 
YCONTRV! 11-6.0 
YCONTRV! 21-0.0 
YCONTRV! 3 ) - YMAXV 
YCONTRV! 41-10.0 

CALL  LGPSPLOT! 1 , ' IXSY’ ,  , YU , 0 , XLABEL , YLABEL ,  ,  ,  .XCONTRU, 

1  YCONTRU,  .UTITLE) 

Setting  up  arrays  for  horizontal  and  vertical  histogram  bar  parameters 

XLOW--200 . 0 
XHI--19S.0 
DO  10  1-1 ,160 
XLOWU1 I ) -XLOW 
XHIUI I ) -XH I 
XLOW-XLOWt-5,0 
XHI-XHI+5.0 
10  CONTINUE 

XLOW— 300.0 
XHI--295 . 0 
DO  20  1-1,120 
XLOWV! I ) -XLOW 


XH I  V (  I  i-XHI 
XLOW-XLOW*  5 . 0 
XHI-XHI*5 . 0 
20  CONTINUE 

r 

C  Plotting  the  horizontal  velocity  histogram 

CALL  LGPSHIST!  1 , XLOWU , XHIU, YU,  160  ,  ,  ,2) 

r 

WRITE1*.*!  'PRESS  <CR>  TO  TERMINATE  PLOT' 

READ! • , ' ( A1 I ’ )  DUMMY 

CALL  LGPS TERM I NATE_PLOT  (  1,1) 

C 

C  Repeating  the  procedure*  to  plot  the  vertical  velocity  hi*togra« 

f 

CALL  LGPSPLOT! 1 , ' IXSY’ ,  . YV , 0 , XLABEL , YLABEL ,  ,  ,  .XCONTRV, 

1  YCONTRV,  , VT I TLE ) 

CALL  LGPSHIST! 1 , XLOWV , XH I V , YV , 1 2  0  ,  ,  ,2) 

C 

WRITE!*,*)  'PRESS  <CR>  TO  TERMINATE  PLOT' 

READ! * , ' ( A 1 ) ' )  DUMMY 

CALL  LGP$TERMINATE_PLOT( 1,1! 

RETURN 

END 


PROGRAM  ORG I 1C 


Written  by  S.  D.  Hedrick  in  July,  1987  to  sequentially  organne  data  files 
to  obtain  a  least  value  of  phase  angle  {third  column  of  data  I  ,  and  associ¬ 
ated  velocities,  at  top  of  file.  This  rearrangement  of  data  was  desired 
to  assist  with  subsequent  manipulation  and  analysis.  In  September,  1987 
this  routine  was  incorporated  as  the  central  program  in  the  processing  and 
analysis  of  unsteady  LDV  velocity  data  generated  by  an  oscillating  airfoil. 
The  program  calls  assisting  subroutines  to  disgard  bad  data  points  and  pro¬ 
vide  wake  profile  plots. 


IMPLICIT  NONE 

REAL  XI  4000)  ,Y( 4000)  ,XPOS, Z  POS , VELO , UVL , WL , UNDIM ( 30 )  , VNDIM!  3  0 ) 
REAL  ZPO!  30 )  , SAV1 , SAV2 , MENU , MENV , PREJ , HEANU ( 18 )  , MEANV (18) 

REAL  MPHAVG ( 18 ) 

INTEGER  214000) , NSAMPLES , PHASA , I , J , X , N , H , SAVE , REJ , STDY 
CHARACTER* 1 0  INPILE 
CHARACTER* 1  ANS 
STDY-0 

WRITE! *,*)' IS  ANALYSIS  POR  NON-OSCILLATING  AIRFOIL?  ( Y/N ) ' 

READ! • , ' ( A 1 ) ’ )  ANS 

IP  (ANS  . EQ-  'Y'  .OR.  ANS  . EQ .  'y')  STDY-1 

WRITE! *,*)' ENTER  NUMBER  OP  DATA  ACQUISITION  POINTS  IN  PROFILE', 
1  <30  maximum  )  ' 

READ!*.*)  N 

WRITE! *,*)' ENTER  PHASE  ANGLE  DESIRED  POR  ANALYSIS  (degrees)' 
READ!*,*)  PHASA 
DO  70  K-l.N 

WRITE! *,*)' ENTER  NAME  OP  INPUT  DATA  FILE  (LDVXXX.DAT)' 

READ! * , ' ( A10 ) ' I  INFILE 

WRITE! *,*)  'ENTER  NUMBER  OP  SAMPLES  IN  DATA  FILE’ 

READ!*,*)  NSAMPLES 

OPEN! UNIT-11 , FILE- INFILE, STATUS- 'OLD' ) 

DO  10  1-1, NSAMPLES 

READ (11,100)  X(I),Y(I),2(I) 

10  CONTINUE 

READ! 11,200)  XPOS.ZPOS 
READ (11,300)  VELO 
IF  (STDY  .EQ.  1)  THEN 
CLOSE! UNIT-11  ! 

CALL  STEADY! NSAMPLES, INFILE, MENU, MENV) 

UNDIM l K ) -MENU /VELO 
VNDIMI  K  f  -MENV /VELO 
GO  TO  80 
END  IF 

M-NSAMPLES-1 
DO  20  J-l, NSAMPLES 
DO  30  I • 1 ,  H 

IF  (Kiel)  .LT.  Z  (  I  )  )  THEN 
SAVE-Z! I ) 

SAVl-X'  I ) 

SAV2-Y!  I ) 

Z( I >-Z( 1*1) 

X!  I ) —  X  <  1*1  ) 

Y! I >-Y< 1*1) 

Z(  I  +  l (-SAVE 
X< I+l l-SAVl 
Y ( I +1 ) -SAV2 
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■  LT-  W-E.  j:-  V-  u  -  t,s.  .1  --  ■ 


£ 


END  ir 

30  CONTINUE 

20  CONTINUE 

CLOSE! UNIT-11 ) 

OPEN! UN IT- 15 , FILE- INFILE , STATUS- 'OLD' ) 

DO  40  I • 1 , NSAMPLES 

WRITE! IS , 100 )  X!  I  )  ,  Y!  I  )  ,  Z(  I  ) 

40  CONTINUE 

CLOSE ( UN IT- 1 5 ) 

WRITE! •,*)' DO  YOU  DESIRE  A  VELOCITY/PHASE  PLOT  OF 
1  '  DATA?  ( Y/N-'DEF-Y ) ’ 

READ! * , ' ( A1 ) ' )  ANS 

IF  (ANS  . EQ .  'N'  .OR.  ANS  . EQ .  'n')  GO  TO  SO 
CALL  PT_PLOT( NSAMPLES , INFILE ) 

SO  CALL  BDPTSI NSAMPLES , INFILE, MENU.MENV, REJ , PREJ , MEA 

1  MPHAVG) 

CALL  AVG! NSAMPLES , INFILE) 

WRITE! *,*)' DO  YOU  DESIRE  INTERPOLATION  OF  THE  DAT. 
READ! • , ' ( A1 ) ' )  ANS 

IF  (ANS  .EQ-  'N'  .OR.  ANS  . EQ .  'n')  GO  TO  55 
CALL  INTERPI NSAMPLES , INTJLE ) 

5S  CALL  CVFI T  (  NSAHPLES  ,  INFILE,  PHASA  ,  UVL  ,  WL  ,  MENU  ,  HEN’. 
1  MEANU , MEANV , MPHAVG ) 

UNDIHI  K  )  -UVL.'VELO 
VNDIM!  K  l-WL/VELO 
80  ZPO(R)-ZPOS 

?0  CONTINUE 

CALL  PROPLOT ( UNDIM , VND I M , ZPO , XPOS , N ) 

Format  statements 

100  FORMAT! SX, F9 . 4 , IX, F9 . 4 . 8X, 14  ) 

200  FORMAT (//,7X,F6.2,7X,F6.2) 

300  FORMAT! 24X, FS . 1 ) 

STOP 

END 


SUBROUTINE  BDPTS ( N SAMPLES , INPILE. MENU, MENV. REJ , PREJ , 

1  MEANU , MEANV .MPHAVG 1 

REAL  X ( 4  000  I  . Y( 4000  )  , MEANU <  18  I  . MEANV!  18  1. STDVU!  18 1 
REAL  STOW  (  18  )  ,  MPHAVG  i  18  1  ,  M I NU  ,  M  I  NV  ,  MAX'J  ,  MAXV  ,  MENU  ,  MENV 
INTEGER  2(4000), PHAS EA , PHAS E B , SUM , COUNT , NSAMPLES , NSAMP . R EJ 
R EAL *  4  XCONTR ( 4 )  , YCONTR I  4  1  , PREJ 
CHARACTER  *  1  ANS , DUMMY , DUMHI 
CHARACTER  *10  INFILE 
CHARACTER  *  1 1  X LABEL 
CHARACTER*21  YLABEU 
CHARACTER* 19  YLABEV 
CHARACTER  *  1 8  TITLE 

OPEN  (  UNIT-1  S , FILE- INFILE , STATUS- ' OLD'  ) 

NSAMP-NSAMPLES 
DO  10  I - 1 , NSAMPLES 

READ! IS, 100)  X(I),Y(I),2(I) 

10  CONTINUE 

100  FORMAT! SX , F9 . 4 , IX , F9 . 4 , 8X, 14 ) 

PHASEA-0 
PHASEB-200 
DO  20  1-1,18 
SUM-0 
SUMU-0 . 

SUMV-0 . 

PHAVG-0 . 

DO  30  J-l, NSAMPLES 

IF  <2(J)  -LE.  PHASEB  .AND.  2(J)  . GE .  PHASEA)  THEN 

SUMU-SUMU-X ( J ! 

SUMV-SUMV- Y  (  J ) 

PHAVG-PHAVG-  (  FLOAT ( 2 ( J )  )/10 .  ) 

SUM-SUM- 1 

END  IF 
30  CONTINUE 

IF  (SUM  .EO.  0)  THEN 

MEANU ( I ) -0 . 

MEANV ( I ) -0  . 

MPHAVG ( I )-( FLOAT ( PHASEA )  + 1 00 . )/10. 

VARU-0 . 

VARV-0 . 

GO  TO  40 
END  IF 

MEANU ( I )-SUMU /FLOAT! SUM) 

MEANV ( I ) -SUMV/ FLOAT! SUM) 

MPHAVG! I )-PHAVG/FLOAT( SUM) 

VARU-0  . 

VARV-0 . 

DO  40  J-l, NSAMPLES 

IF  ( 2 ( J )  . LE .  PHASEB  .AND.  2(J)  .GE.  PHASEA)  THEN 
A1 -X ( J ) -MEANU ( I ) 

A2- Y( J ) -MEANV { I ) 

B l -Al *  *  2 /FLOAT ( SUM ) 

B2-A2  *  *  2 /FLOAT ( SUM ) 

VARU-VARU+Bl 

VARV-VARV+B2 

END  IF 
40  CONTINUE 

STDVU( I )-SQRT( VARU ) 

STDWf  I  )-SQRT<  VARV) 

PHAS EA- PHAS EA- 2 0 0 
PHASES- PHASEB- 2 00 
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20  CONTINUE 

WRITE! • , • I ' DO  YOU  DESIRE  A  STANDARD  DEVIATION  PLOT?  ( Y/N  -  DEF-Y)' 
READ!  •  ,  200 )  ANS 
200  PORRAT ( A1 ) 

IF  (ANS  . EO ■  'N'  .OR.  ANS  . EQ .  'n')  THEN 

GO  TO  50 
END  IT 

Setting  up  the  coordinate  axis  lor  the  plot 

XLABEL- ' PHASE  ANGLE' 

YLABEU-' HORIZONTAL  VELOCITIES' 

YLABEV- ' VERTICAL  VELOCITIES' 

TITLE- ' STANDARD  DEVIATION' 

XCONTR ( 11-6.0 
XCONTR ( 21-0.0 
XCONTR! 3  t - 360 . 0 
XCONTR! 41-45.0 
YCONTR I  11-6.0 
YCONTR ( 2  >--200 . 0 
YCONTR! 3 1 -600 . 0 
YCONTR! 4 1-100 .0 

CALL  LGPSPLOT! 1 .' EXSY ' .HPHAVG, HEANU, 0 .XLABEL, YLABEU ,  ,  ,  .XCONTR, 

1  YCONTR,  .TITLE) 

Plotting  the  horizontal  velocity  data  points  on  the  axis 

CALL  LG PS  PO I NT 1  1 , HPHAVG , KEANU ,  .03,4,16,  ,  ,  ) 

Plotting  the  horizontal  standard  deviation 

CALL  LGPSSTNDEV! 1 , HPHAVG. HEANU, STDVU, 18  ,  1 

45  WRITE! *,•)  'PRESS  <CR>  TO  TERHINATE  PLOT' 

READ! *,200)  DUMMY 
CALL  LGPSTERMINATE_PLOT! 1,1) 

Plotting  the  vertical  standard  deviation 

CALL  LGPSPLOT! 1 , ' EXSY' .HPHAVG, MEANV.0, XLABEL, YLABEV,  ,  ,  .XCONTR, 

1  YCONTR,  .TITLE) 

CALL  LGPS  PO I  NT ! 1 , HPHAVG , HEANV ,  .03,4,16,  ,  ,  ) 

CALL  LGPSSTNDEV!  1  .HPHAVG. MEANV.  STDW,  18,  1 

47  WRITE! *,*>  'PRESS  <CR>  TO  TERHINATE  PLOT' 

READ! •, 200)  DUHHI 
CALL  LGPSTERHINATE_PLOT( 1,1) 

C 

50  COUNT-1 
PHASEA-0 
PHASEB-200 
DO  60  1-1,18 

HINU-HEANU! I ! -2 . 0» STDVU! I ) 

HINV-MEANV!  I  )  -  2 .  0*STDW(  I  ) 

RAXU-HEANU ( I ) ♦ 2 . 0 • STDVU ( I ) 

MAX V- HEANV (  I  )*2 . 0»STDW!  I  ) 

DO  70  J- 1 , NSAHPLES 

IP  (Z(J)  .LE.  PHASES  .AND.  ZU>  . GE .  PHASEA)  THEN 
IP  (X(J)  .  LT .  HINU  .OR.  XU)  .GT.  RAXU )  THEN 
COUNT-COUNT* 1 
IP  ( J  .EQ.  1 )  GO  TO  70 


80 


85 

70 

60 


90 


95 


97 


98 

96 


00  80  X-O.J-2 
X( J-X  > - X ( J-X-l  I 
Yl J-X  I  - Y  t J-X-l ) 

Z!  J-X l-Zl J-R-l  ) 

CONTINUE 
GO  TO  70 
END  IF 

IF  (  Y  !  J  )  .LT.  MINV  .OR.  Y(J)  . GT .  MAXV )  THEN 
COUNT-COUNT* 1 
I F  !  J  .  EQ .  II  GO  TO  7  0 
DO  85  X-O.J-2 
X: J-X  )-X( J-X-l I 
Y( J-X  > - Y ( J-X-l ) 

Z( J-X )-Z( J-K-l ) 

CONTINUE 
END  If 
END  If 
CONTINUE 

PHASER- PHASER* 2 00 
P HAS EB- PHASES*  200 
CONTINUE 
CLOSE! UNIT-1  5  I 

OPEN! UNIT-12 , f I LE- INTI LE, STATUS- ' OLD' ) 

DO  90  I  -  COUNT , NS AH PL ES 

WRITE! 12 , 100 )  X!I),Y(I),Z(I) 

CONTINUE 

NSAMPLES-NSAMPLES* 1 -COUNT 
CLOSE ( UN  I T- 1 2  ) 

OPEN (UN IT-13, FILE-INFILE, STATUS- 'OLD' ) 

SUU-0 . 

SUV-0 . 

DO  95  1-1 , NSAMPLES 

READ!  13.100)  X!  I  )  , Y( I ) ,Z(  I) 

SUU-SUU*X(  I  ) 

SUV-SUV*Y(  I  I 
CONTINUE 

MENU-SUU/FLOAT! NSAMPLES ) 

MENV-SUV/ FLOAT! NSAMPLES) 

CLOSE ( UN I T- 1 3 ) 

REJ-NSAMP-NSAMPLES 

PREJ-! 1 .-<  FLOAT) NSAMPLES  t/FLOAT(NSAMP) 1) *100 . 

PHASEA-0 
PHASEB-200 
DO  96  1-1.18 
SUM-0 
SUHU-0 . 

SUHV-0 . 

PHAVG-0 . 

DO  97  J-l .NSAMPLES 

IF  (Z(JI  .LE.  PHASEB  .AND.  Z(J)  . GE .  PHASER)  THEN 
SUHU-SUHU+X) J ) 

SUHV-SUHV*Y! J 1 

PHAVG-PHAVG*) rLOAT) Z) J! )/10.  ) 

SUM-SUM* 1 
END  IF 
CONTINUE 

IF  (SUM  .EQ.  01  THEN 
MEANU! I  1-0. 

MEANV ( I )-0. 

MPHAVG ( I )-( FLOAT) PHASER  1*100.  )/10. 

GO  TO  98 
END  IP 

MEANU) I  I -SUHU/FLOAT) SUM) 

MEANV) I ) — SUM V/ FLOAT ( SUM) 

MPHAVG) I  I - PHAVG/FLOAT l SUM ) 

PKAS  EA- PHAS  EA*  2  0  0 
PHASEB- PHAS EB* 2 00 
CONTINUE 
RETURN 
END 


63 


.  . 


-VlV-JL 


CmJsLrnLi. 


SUBROUTINE  AVG I NSANPLES . INFILE) 

REAL  X( 4000  )  .  Y  4000  i 

INTEGER  Z ( 4000 ) .COUNT , SUN , NSANPLES 
CHARACTER  *10  INFILE 
COUNT-1 
N-NSAMPLES- 1 

OPEN. UNIT-1 FILE. INFILE , STATUS- 'OLD' ) 

DC  10  I  - 1 . NSANPLES 

READ ( 12.100)  X<  I ) , Y( I) ,Z( I  ) 

10  CONTINUE 

100  FORNATi SX, F9 . 4 , IX, T9 . 4 , 8X, 14 > 

DO  20  I - 1 , N 

IF  (Z(I.l)  .EQ.  2(1)  .AND.  2(1.2)  . NE .  2(1))  THEN 

COUNT- COUNT. 1 
X(I.l)-(X(I*l).X(I))/2.0 
Y< 1.1 )-( Y( I.l ).Y( I ) )/2 . 0 
DO  30  J-0.I-2 
X(  I-J )-X(  I-J-l  I 
Y ( I-J)-Y( 1-J-l) 

Z( I-J )-Z(  I-J-l  ) 

30  CONTINUE 

END  ir 

IF  (2(1.2)  .EQ.  Z(I))  THEN 
SUM-1 

DO  25  J- 2 , NSANPLES- I 

IF  (Z(  I+J)  .EQ.  Z  l  I  )  )  THEN 
SUN-SUN. 1 
END  ir 
25  CONTINUE 

DO  27  J-l.SUH 

X(  I  +  SUN ) - X ( I *SUN ) *X(  I+SUN-J ) 

Y ( I ♦ SUN ) - Y ( I +SUN I + Y(  I.SUN-J) 

27  CONTINUE 

X(  I .SUM ) -X (  I ♦ SUN  I / ( SUN* 1  ) 

Y (  I *SUN l-Y( 1 *  SUM ) / ( SUN. 1 ) 

DO  28  J. 1,1-1 

X< I.SUN-J  )-X( I-J ) 

Y (  I.SUN-J  )-Y( I-J ) 

Z( I .SUM- J ) -Z ( I-J) 

28  CONTINUE 
COUNT- COUNT* SUN 

END  IF 
20  CONTINUE 

CLOSE( UNIT-12) 

OPEN! UNIT- 1 4 , F I LE- I NT I LE , STATUS- ' OLD' ) 

DO  40  I -COUNT, NSANPLES 

WRITE! 14,100)  X(I),Y(I),Z(I) 

40  CONTINUE 

NSANPLES-NSANPLES+1 -COUNT 
CLOSE! UNIT-14  ) 

RETURN 

END 


) 


SUBROUTINE  INTERPI NSAMPLES, INFILE) 

IMPLICIT  NONE 
CHARACTER  *10  INFILE 
REAL  XI 3600 ) , Yl 3600 ) 

INTEGER  1(3600)  , NSAMPLES , I , J , COUNT, NS A, FAC 

CHARACTER • 1  ANS 

OPEN) UNIT-11 , FILE- INFILE , STATUS- 'OLD' ) 

DO  10  1-1, NSAMPLES 

READ! 11,100)  X(I),Y(I),Z(I) 

CONTINUE 

FORMAT) 5X,F9.4,1X,F9.4,8X,I4) 

COUNT- 0 
NSA-NSAMPLES 
DO  20  1-1,3600 

IF  (COUNT  ,GE.  NSA!  THEN 
Z  (  I  )  -  I  -  1 

X(  I  )-(  (X( 1  ) —  X (  I  - 1 )  )/(2( 1 >-Z(  1-1  )  )  ) *  ( Z(  I  )-Z(  1-1  )  )  *X  (  1-1  ) 

Y(  I  )-(  (Y(  1  )-Y(  I  —  1 )  )/(Zll )  —  2 (  1-1 )  )  ) * ( 2 ( I )  —  2 (  1-1  )  )  ♦  Y  (  1-1  ) 

NSAMPLES-NSAMPLES* 1 
GO  TO  20 
END  IF 

IF  (2(1)  . EQ .  I - 1 )  COUNT-COUNT* 1 

IF  (Z(I)  .NE.  I - 1 )  THEN 
NSAMPLES-NSAMPLES*! 

DO  30  J- 0 , NSAMPLES- I  - 1 

X( NSAMPLES- J * -X ( NSAMPLES -J-l ) 

Y ( NSAMPLES- J )- Y ( NSAMPLES- J- 1 ) 

2 ( NSAMPLES- J ) -Z I NSAMPLES- J - 1 ) 

CONTINUE 
Z (  I  )-I-l 

IF  (I  . EQ .  1)  THEN 

X(t)-!(X(t*l  )-X(  NSAMPLES  )  )/(Z(I*l)-Z(  NSAMPLES  )  )  )  *  (  2  (  2  >- 
2 ( NSAMPLES I  I  *X( NSAMPLES > 

Y ( I ) - ( ! Y ( I ♦ 1 ) - Y ( NSAMPLES ) */(Z( 1*1 )-2( NSAMPLES ) ) I • ( 2! I ) - 
2  (  NSAMPLES )  > *Y( NSAMPLES ) 

GO  TO  20 
END  IF 

X!Ii-((X(I*l)-X(I-l))/(Z(I*l)-Z(I-l)))*(Z(I)-Z(I-l))*X(I-l) 
Y(  I  1 •  (  1  Y(  1*1 ) —  Y (  1-1  )  )/(Z(  1*1  )-Z(  1-1  )  )  )*( 2 (  I  )-Z(  1-1  )  )  ♦  Y  (  1-1  ) 
END  IF 
CONTINUE 
CLOSE! UNIT-11 ) 

WRITE) *,*)' DO  YOU  DESIRE  TO  REDUCE  THE  NUMBER  OF  INTERPOLATED', 
1  '  DATA  POINTS  BY  AN  INTEGER  FACTOR  OF  3600?  IY/N)' 

READ' * , ' I A1 ) ' )  ANS 

OPEN! UNIT-12 , r I LE- INFILE, STATUS- 'OLD' ) 

IF  (ANS  .EQ.  'Y'  .OR.  ANS  . EQ .  'y')  THEN 
WRITE! •,*)' ENTER  INTEGER  FACTOR' 

READ)*,*)  FAC 
FAC-20 

DO  34  I-FAC , NSAMPLES , TAC 

WRITE!  12,100)  X( I ( , Y<  I )  , Z(  I  ) 

CONTINUE 

NSAMPLES-3600/FAC 
GO  TO  40 
END  IF 

DO  40  I- I, NSAMPLES 

WRITE) 12 , 100 )  X(I),Y(I),ZII) 

CONTINUE 
CLOSE! UNIT-12  ) 

RETURN 

END 
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SUBROUTINE  FOURI ER (  F , H , A , 8  ,  N  ) 

DIMENSION  F(0:M), FA (0:300 ),FBI0:300l 
DATA  PI/3  .  14159265/ 

H»2*PI/M 
TYPE*,  H 
J-0 

DO  2  I-0,M 
04 

IF  (I  .  EQ .  J  !  02 
IT  (I  .EQ.  Jl  002 
FA( I  )  —  F (  I  > *COS I N • H *  I  ) 

FB( I >-F{ I ) *SIN( N*H* I ) 

IF  (I  . EQ .  Ml  T3-FA( I ) 

IF  (I  . EQ.  M)  TT3-FBI  I ) 

IF  (I  -EQ.  M-l)  T2-FAII) 

IF  (I  .EQ.  M-l)  TT2-FBII) 

IF  (I  .EQ.  H-2)  Tl-FA(I) 

IF  (I  .EQ.  M-2)  TTl-FB(I) 

IF  (I  .EQ.  M-3)  TO-FAI I ) 

IF  (I  .EQ.  M-3)  TTO-FBII) 

IF  (  I  . EQ .  0  I  GO  TO  2 

FA(  I  )-FA(  1-1  )«-C*FA( I ) 

FBI  I  ) - FB (  1-1 ) +C • FB ( I  ) 

CONTINUE 

A-H/3./PIM  rA(M)-T3) 

B-H/3 . /PI  * ( FBI M ) -TT3 ) 

IF  I  J  .  NE .  M-*l  )  GO  TO  3 

A-H/3 . /PI • ( FAI M)-4 . *T3-2 . *T2-4 . *T1-T0 ) 

A- A*  3. *H/8 ,/PI*(T0*3. *T1 ♦ 3 . *T2  *T  3  3 
B-H/3./PI*(FB(M)-4.*TT3-2. «TT2-4 . *TT1-TT0 ) 
B-B*  3 .  *H/8 ./PI • I TT0  +  3 . *TT1  ♦  3 .  *TT2  +  TT 3 ) 
CONTINUE 
RETURN 


r  r  r rr r *  a 


SUBROUTINE  CVTITINSAMPLES, INFILE, PHASA, UVL, WL, MENU,  MENV.  REJ  ,  PREJ  , 
1  MEANU . MEANV . MPHAVG ) 


written  by  S.  D.  Hedrick  in  October,  198"'  to  fit  a  sinusoidal  curve  to  LDV 
velocity  data  aenerated  by  an  oscillating  airfoil.  The  subroutine,  utiliz¬ 
ing  FOURIER,  adapted  from  June  19"’J  DCS  Airfoil  program,  fits  a  single 
cycle  sine  wave  to  tbe  data  and  computes  the  u-  and  v-velocity  components 
at  a  phase  angle  specified  for  the  analysis. 


IMPLICIT  NONE 
CHARACTERS  ANS  .  DUMMY 
CHARACTER  *10  INFILE 

INTEGER  NSAMPLES , I , PHASA ,NSA, REJ ,ARG 

REAL  FU!  0 :  36  00  )  , FV( 0 : 36  0  0 i ,AU, AV , BU, BV. PI , X(  36  00  )  ,  Y(  36  00  ) 

REAL  UVL , WL . MENU , MENV , PREJ , MEANU ( 18 ) , MEANV' 18 ) , MPHAVG1 18 ) 

DATA  PI  '3.  1  4  1  592  7/ 

WRITE' * , * ) 'DESIRE  A  VELOCITY/PHASE  PLOT  OF  THE  REFINED  DATA?' 
WRITE! *,•)'( Y/N--DEF-Y) ' 

READ!  • ,  '  t  A1  t '  )  ANS 

IF  (ANS  .EQ-  'N'  .OR.  ANS  . EQ .  'n't  GO  TO  7 
CALL  RPLOT I NSAMPLES. INFILE, REJ , PREJ 1 
OPENi UNIT-11 , FILE-INF  I LE , STATUS- ' OLD'  ) 

DO  5  I-i, NSAMPLES 

READ!  11  ,  100  I  X (  I  I  , Y (  I  I 
CONTINUE 

FORMAT'  5X . F9 . 4 , IX , F9 . 4  ) 

CLOSE' UNIT-11 i 
DO  10  1-1. NSAMPLES 
FU  •  1  -  1  '  -  X  I  I  ) 

FVi  I  - 1  ) - Y (  I  ) 

CONTINUE 

NSA-NSAMPLES-1 

CALL  FOURIER' FU.NSA.AU, BU,1 ) 

CALL  FOURIER! FV,NSA,AV, BV, 1 ) 

DC  20  I-0.NSA 

FUI  I *1  >-M£NU*AU*COS (  ( 2 . *  PI /FLOAT' NSAMPLES )  I • FLOAT!  I  )  ) ♦ 

1  BU*  SIN'  ( 2 . *PI /FLOAT! NSAMPLES i ' » FLOAT' 1 ) ) 

FV (  I ♦ 1  ) -MENV* AV  *  COS (  ( 2 . *PI /FLOAT! NSAMPLES I  1 ‘FLOAT!  I  )  ) ♦ 

1  BVSINi  (  2  .  *PI/rLOAT(  NSAMPLES  I  1  *  FLOAT'  I  )  ) 

CONTINUE 

WRITE' *,* l ’ PLOTS  WILL  NOW  BE  PROVIDED  TO  COMPARE  SINUSOID  WITH' 
1  '  VELOCITY  DATA.  PRESS  <CR>  TO  CONTINUE.' 

READ' • , ' ( A1 I ' )  DUMMY 

CALL  PLOT! NSAMPLES, I NF I LE , PU , FV , MENU , MENV , MEANU , MEANV , 

1  HPHAVG) 

ARG-NSAMPLES* ( FLOAT! PHASA1/360. ) 

WRITE! • . • I '  ARG  -' , ARC 
UVL- FU ' ARC  I 
WL-FVI  ARG  ) 

RETURN 

END 


SUBROUTINE  PLOT! NSAMPLES , I NFI  LE , FU , EV , MENU ,  MENV  , MEANU, 
1  MEANV , MPHAVG > 


Written  by  S.  0.  Hedrick  in  October,  1987  to  plot  a  comparison  of  velocity 
data  and  a  sinusoidal  curve  fit. 


IMPLICIT  NONE 

INTEGER  NSAMPLES , I , BOUND, STOPU , STOPV 

REAL* 4  XVAL <  4000  1  , YVALUi  4  000 1  , YVALV (  4  0  00 )  , XCONTR!  4 )  , YCONTU!  4  ) 
REAL *4  YCONTV! 4 ) , FU! NSAMPLES ) , FV( NSAMPLES ) 

REAL  *4  YHINU,  YMAXU  .  YHINV,  YMAXV  ,  MENU  .  MENV 
REAL  MEANUt 18 ) , NEANV ( 18 1 ,MPHAVG( 18 ) 

CHARACTER  *  1 0  INFILE 
CHARACTER* 1 1  X LABEL 
CHARACTER*15  YLABEL 
CHARACTER  *  2 1  UTITLE 
CHARACTER  *  1 9  VTITLE 
CHARACTER  *  1  DUMMY 

OPEN! UNIT-11 . FILE- INFILE, STATUS- 'OLD' ) 

DO  10  I - 1 , NSAMPLES 

READ! 11,100)  YVALU ( I ) , YVALV ! I ) ,XVAL( I ) 

XVAL!  I  I -XVAL ( I  1/10. 

10  CONTINUE 

100  TORMAT! 5X,F9.4,1X,F9.4,8X,F4.0) 

XLABEL- ' PHASE  ANGLE' 

YLABEL- ' VELOCITY  (ft/S)' 

UTITLE- ’ HORIZONTAL  VELOCITIES' 

VTITLE- ' VERTICAL  VELOCITIES' 

YMINU-MENU-1S0 . 

YMAXU- MENU* 150 . 

YM I NV-HENV- ISO . 

YMAXV-MENV-150 . 

XCONTRi  11-6.0 
XCONTR! 21-0.0 
XCONTR! 3 1-360 . 0 
XCONTR! 41-45.0 
YCONTU!  1  ) - 6  -  0 
YCONTU1 2 l-YHINU 
YCONTU! 3 ) -YMAXU 
YCONTU'  4  i-50 . 0 
Y  CONTV i 11-6.0 
YCONTV  2 i-YMINV 
YCONTV 1 3 i -YMAXV 
YCONTV  I  4'-50.0 

CALL  LGP$ PLOT ( 1 .  ' EXS Y '  , XVAL , FU , NSAMPLES , XLABEL , YLABEL ,  ,  ,  , 

1  XCONTR, YCONTU,  .UTITLE) 

CALL  LGPSPOINT!  1 .MPHAVG, MEANU,  . 03  ,  3  ,  18  ,  ,  ,  ) 

WRITE'  *,*  I  ' PRESS  <CR >  TO  TERMINATE  PLOT' 

READ! *,200)  DUMMY 

CALL  LG PS TERM I NATE  PLOT! 1.1) 

CALL  LGPSPLOT! 1, 'EXSY' , XVAL, FV, NSAMPLES, XLABEL, YLABEL,  ,  ,  , 

1  XCONTR, YCONTV,  .VTITLE) 

CALL  LGPSPOINT! 1 , MPHAVG , NEANV , . 03 , 3 , 18 ,  ,  ,  > 

WRITE!*,*)  'Pteaa  <CR>  to  terminate  plot' 

READ! *.200)  DUMMY 
200  FORMAT! Al ) 

CALL  LGPSTERMIMATE  PLOT! 1,1) 

CLOSE ( UNIT-1 1 ) 

RETURN 

END 
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SUBROUTINE  PRO PLOT! UNDIH , VND  IN , 2  PO , XPOS , N I 
IMPLICIT  NONE 

R  EAL • 4  UNDI M '  30  , VNDI M '  3C ' . 2 PO '  30 ' . XCONTU  4  '.XCONTV I  4 
REAL*  4  YCONTR!  4  ,  , UNDI !  3000  ,VNOI  3 C C 0  i  , 2 P  •  3 0 0 0 i 

REAL  XPOS 
INTEGER  N.NBR 


CHARACTER  *29 

XLABEU . XLABEV 

CHARACTER • 3  2 

YLABEL 

:haractep* : ? 

UTITLE 

CHARACTER *2S 

VTITLE 

CHARACTER • 1 

D'JHHY . ANS 

CHARACTER'S 

XPO 

chapattep  •  4 

XXPO 

XLABEU- • VELv 

r  :  TY  RATIO  *J  T 

XLABEV. ■ VELOCI TY  RATIO  ‘V  UINTIN 
YLA8EL. ' PROBE  POSITION  IN.  FROM  CENTER'1 
UTITLE-  HORIZONTAL  VELOCITY  PROP  I LE  * 

VTITLE- ’ 7ERTI CAL  VELOCITY  PROrtLE' 

XPOS-XPOS  3  0 

WRITE  UNI T-XXPC . FMT- ' i T4 . 1 ' ' I  XPOS 
XPO-  X  C  -■  XXPO 

XCONTU  I '-6.0 
XCONTU  ; > - C  0 
XCONTU  3'-:. 2 
XCONTU  4 i -0 . 2 
XCONTV ■ 1 . - e . 0 
XCONTV  2  ■  --C . 7 
XCONTV  3  '  -  0  .  7 
XCONTV.  4  -C  .  2 
YCONTP  1-6.0 
YCONTR  2  -  -  2  C 
YCONTP  3 '-2.0 
YCONTR  4 
NBP-10*  N -  1  *1 

CALL  LOPS  PLOT  1  .  ' EXS Y '  . UNO  I M . 2PO , 0 , XL AS EU . YLABEL , 

1  ,  .  .XCONTU. YCONTR,  , UTITLE) 

CALL  LOPS  POINT  1 , UNO I M . 2 PO ,  . 0 3 , 3 , N ,  ,  ,  ) 

CALL  LGPSPUT  TEXT' 1 . . 4 , 5 . 5 . XPO,  ,  ! 

WP.'TE  *.*  PRESS  <  CR  >  TO  TERMINATE  PLOT' 

READ  • . ■  A1  ' •  DUMMY 

CALL  LOPS  TERM  I  NATE  PLOT'  1,1) 

CALL  LGPSPLOT  2 .' EXS Y '. VND I H , 2PO . 0 , XLABEV , YLABEL , 

1  .  .  , XCONTV, YCONTR,  .VTITLE) 

CALL  LOPS  PCI  NT  1  . VNO I M . 2 PO .  . 0 3 , 3 , N ,  ,  ,  t 

CALL  LGPSPUT  TEXT  1 . . 4 . S . 5 . XPO .  .  I 

WRITE  *.*  PRESS  <CR>  TC  TERMINATE  PLOT' 

READ' • , ' • AI  '  DUMMY 

CALL  LGPSTERMINATE  PLOT! 1 . 1 ) 

WRITE1 *,*  ’DO  YOU  DESIRE  A  SPLINE  INTERPOLATION  CURVE  T 
1  '  THE  DATA  POINTS?  ( Y  N ' ' 

READ' • , ' I Al > '  ANS 

IP  IANS  . EQ .  'N'  .OR.  ANS  . EQ .  'n'  >  GO  TC  10 

CALL  LGPSPLOT! 1. 'EXSY' .UNDIM.2PO.O . X LABEL . YLABEL . 

1  .  ,  , XCONTU , YCONTR ,  .UTITLE' 

CALL  LGPSPOINTl  1 . UNDIM. 2PO ,  . 0 3  .  3  ,  N ,  .  , 

CALL  LGPSSPLINE  QHC ( UNO IM.2PO.N.10, UND I  , 2  P .  ) 

CALL  LGPS  PLOTM ( I ,’EXSY'.UNDI.2P,NBR,I 
CALL  LGPSPUT  TEXT ( 1 , . 4 . 5 . $ . XPO .  .  1 

WRITE! *,* I ' PRESS  <CR>  TO  TERMINATE  PLOT' 

READ! • , ' ( Al I ' )  DUMMY 

CALL  LGPSTERM I  NATE _ PLOT '  1  ,  1  I 

CALL  LGPSPLOT! 1 , 'EXSY' . VND I M , 2 PO . C , X LAB EV . YLABEL , 

1  ,  ,  .XCONTV. YCONTR.  .VTITLE) 

CALL  LGPSPOINTl 1 , VNDIM , IPO, . 03 . 3 , N,  ,  ,  ) 

CALL  LGPSSPLINE  QHC ( VND I M . Z PO , N ,  1 0 , VND I  . 2 P  ,  ) 

CALL  LGPS PLOTM! I ,  'EXSY'  , VND I  , ZP , NBR ,  ) 

CALL  LGPSPUT  TEXT  I  1 ,  . 4 . S . 5 . X PO ,  ,  ) 

WRITE! *,* I ' PRESS  <CR>  TO  TERMINATE  PLOT' 

READ!  • ,  '  ( Al  )  '  )  DUMMY 

CALL  LGPS TERM  I NATE_ PLOT ( 1,1) 

10  RETURN 
END 
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